MICRO-CARO 


i  Reproduce  by 


CONTROL  ONLY 


ARLINGTON  HALL  STATION;  ARLINGTON  12  VIRGINIA 


ll-- 

Ip  ' 

i ' 

n 

1  r^B 

1 

H'  r 

(  '¥ 

-  ? 

\  F 

P  I 

r  «  !■ 

H 

'1 

1'  *  ^ 

1 

■ ' 

/  1 

1  j| 

1 

r 

1  _  i_  J 

"NOTICE:  When  Government  or  other  drawing*,  specification*  or 
other  data  are  ust^d  for  any  purpose  other  than  in  connection  wit 
a  definitely  related  Government  procurement  operation,  the  U.S. 
Government  thereby  incur*  no  responsibility,  nor  any  obligation 
whatsoever:  and  the  fact  that  the  Government  may  have  formulate 
furnished,  or  in  any  way  supplied  the  said  drawings,  specification* 
or  other  data  is  not  to  be  regarded  by  implication  or  otherwise  a 
in  any  manner  licensing  the  holder  or  any  other  person  or  corpor. 
tion,  or  conveying  any  right*  or  permission  to  manufacture,  use  or 
sell  any  patented  invention  that  may  in  any  way  be  related  theret' 


J 

J 


r 


UNIVERSITY  CALIFORNIA 


INSTITUTE  OF  EN 
BERKELE 


ERING  RESEARCH 
LIFORNIA 


*• 

Ml 


WATER  WAVES.  IV. 


7 

C. 

o 

I 

X 

>  *•" 
r  > 


With  partial  support  from 
Contract  Nuabar  N-onr-222(30) 


SERIES  NO. 

ISSUE  NO 
DATE 

! 


82 

II 

Wovambtr  1959 


f\ll  COPY 


WATER  WAVES.  IV 


by 

E.«  V.  l4iitone 


with  partial  support 
from 

Contract  Nwbar  Il-onr-222(30) 


Unlvarslty  of  California 
Institute  of  Enginatring  Raaaarch 
Barkalay  4,  California 


Fovsn^ar,  1959 


Xliis  report  coastltutcs  the  fourth  pert  of 
an  article  on  Water  Waves  being  prepared  for  the 
new  edition  of  the  Encyclopedia  of  Physics  (Handbuch 
der  Physik)  published  by  Springer  Verlag. 


Shallow  watar  . . . . 

23.  The  fvmdamantal  aquations  for  tha  first  approjcimatlon  oE-2 

29.  Tha  llnaarlaad  shallow-watar  theory . a  .  .  E-4 


290C.  Unaarlaad  shallow-water  theory  applied  to 
two-diaensional  steady  flow  .  . . 


30.  Honllnaar  shallow-water  tlieory 


2-17 


30 dC.  One-dlaenslcnal  non-steady,  first-order, 

shallow-water  theory . . . E-17 

30^.  Two- dinanslonal,  steady,  supercritical  flow  by 
•  the  first-order  shallow-water  theory . E-24 

30 yt  One-dimensional,  steady,  open  channel  hydraulics 
and  the  hydraulic  jump . . 

31.  Higher  order  theories  and  the  solitary  cind  cnoidal 

waves  . . . . 


31(/.  The  first  and  second  approximations  to  the 
cnoidal  find  solitary  waves . . 


E-58 


The  limiting  height  and  velocity  of  propagation 

of  cnoidal  and  solitary  waves  .  . 

» 


E-1 


2.  Shallow  Watar  Waves 

Unis  chapter  %d.ll  deal  with  special  solutions  based  on  the 
shallow-water  approxlj&atlon»  following  the  method  of  Friedrichs 
[1948]  as  presented  In  section  10/&  *  The  shallow-water  approxima¬ 
tion  for  the  waves  over  a  rigid  bottom  yields  a  set  of  nonlinear 
equations  (cf*  10.32)  Sven  In  the  first  approximation.  If  these 
equations  are  then  linearized  they  result  in  a  hyperbolic  type 
equation  which  reduces  to  the  simple  wave  equation  for  a  flat 
horizontal  bottom.  Consequently  the  solutions  resulting  from  the 
shallow-water  approximation  are  completely  different  in  character 
from  those  resulting  from  the  Influlteslmal-wave  approximation  of 
section  10 OC  and  chapter  D,  which  resulted  In  linear  equations  and 
linear  boundary  conditions.  That  l8»  the  shallow-water  approximation 
leads  to  nonlinear  hyperbolic  type  equations,  whereas  the  Inflnltes- 
imal-wave  approximation  leads  to  a  set  of  linear  equations  satisfying 
the  boundary  conditions  and  having  each  successive  approximation  to 
the  velocity  potential  satisfy  the  simplest  elliptic  equation, 
nsMly  the  Laplace  equation. 

After  the  first-order  shallow-water  approximation  (10.32)  has 
been  applied  to  several  problems,  then  the  method  of  Friedrich::  [1948] 
and  Keller  [1948]  will  be  extencied  to  obtain  the  second  and  third 
approximations  of  the  shallow-water  theory  and  thereby  present,  for 
the  first  time,  the  exact  second  approximation  to  the  cnoldal  wave 
of  Korteweg  and  de  Vries  [1895],  and  the  solitary  wave  of  Bousslnesq 
[1871],  and  Rayleigh  [1876].  These  higher  order  approximations  lead 
directly  to  relations  predicting  the  aaxlmtai  heights  of  cnoldal 
waves  and  solitary  waves. 


I 
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28.  Th«  Fun<UMntal  Equations  for  the  First  Approximation. 

The  shallow-water  expansion  metihod  introduced  by  Friedrichs 
[1943]  is  discussed  in  section  10.  For  this  application  the  ex¬ 
pansion  parameter  S  was  selected  so  that  the  first  approximation 
would  be  identical  to  the  nonlinear  equations  of  the  classical 
shallow-water  theory >  which  is  based  on  the  assumption  of  hydro- 
gtatic  pressure  variation  throughout  and  the  neglect  of  the  variaton 
of  the  horizontal  velocity  components  with  depth,  so  the  complicated 
boundary  value  problem  is  greatly  simplified  to  the  following  non¬ 
linear  equations: 

Ut  +  -  -gy^ 

+  uWjj  +  wWj  -  (28.1) 

t  +  t'i(7+  b)]^  +  [w(7  +  h)]jj.  -  0 
See  Lamb  [1932],  page  254,  or  Stoker  [1957],  page  23. 

Upon  noting  (see  Figure  37)  that  the  rigid  hotter:'  surface  or 
the  undisturbed  water  depth  in  (28.1)  is  designated  by 
y— h(x,z)-b(x,z)  for  the  corresponding  designation  in  (10.32),  we 
see  that  the  (28.1)  set  of  nonlinear  equations  are  identical  to 

(10.32)  which  correspond  to  the  first  approximation  of  Friedrichs 
[1948]  shallow-water  expansion  method  as  discussed  in  sectirm  10, 
thereby  showing  that  the  method  has  physical  Justification.  It  is 
seen  that  the  higher  order  approximations  following  (10.25)  and 

(10.33)  will  require*  that  £  be  sufficiently  small,  consequently, 
as  will  be  shown,  this  expansion  method  is  applicable  only  if  the 
water  depth  and  surface  curvature  are  small  in  comparison  to  the 
equivalent  of  tho  wave  length  of  the  free  surface  profile  disturbance. 


E-3 


IT 

ir 
>» 

f 

a 

Thftrt£ora»  in  som  caaas,  this  shallow  watsr  thaory  is  applicsbls 
t:o  axtrswaly  !«?(•  wstsr  dspths  ss  long  as  the  wavn  length  is 
sufficiently  long,  the  most  ccaiaon  application  being  to  tidal 
wives  which  are  the  oceanic  tides  produced  by  the  gravitational  action 
of  the  sun  and  the  moon,  e*g*  see  Lamb  [1932]* 

The  isatheoutical  Justification  for  this  shallow-water  expansion 
method,  at  least  for  special  cases,  lies  in  the  existence  proof  of 
Friedrichs  and  Hyers  [1954]  for  the  solitary  wave,  and  the  existence 
proof  of  Littman  [1957]  for  the  more  general  cnoidal  waves*  Both 
of  these  proofs  demonstrate  that  this  expansion  method  will  at  least 
give  as3miptotio  discripticms  of  the  exact  solutions  for  these 
particular  problems* 

The  nonlinear  first  approximation  given  by  (28*1)  is  considerably 
simplified  if  the  rigid  bottom  surface  h(x,s)  is  flat  and  horisontal, 

ms  my  be  seen  by  letting  h  •  constant  so  that  (28*1)  may  be 

written  as  +  udjj  +  wu^  -  -g(^  +  h)^^ 

Wt  +  ler*  +  -  -g('J^  +  h)^  (28*2) 

+  h)t  +  lu(^  +  h)]j^  +  (w(«^  +  h)]^  -  0 

which  is  identical  to  the  well  known  two- dimensions  I  gas  dynamics  ) 

equation,  e*g,  see  Lamb  [1932],  if  we  write 

^  (x,z,t)  -  (x;Z,t)  +  h]  (28*3) 

yp  -  « 

pt.  o 

Since  the  isentropic  gas  relationship  is  p  *  (constant.]^"  ), 
therefore  the  first  order  nonlinear  .*ihallow-vater  approximation 
for  I  flat  horizontal  bottom  is  Identical  to  thft  isentropic  two- 
dlmensioiial  gas  flow  having  a  specific  heat  ratio  of  2*  This 
Is  the  basis  of  the  so  called  hydraulic  analogy  which  has  been 


constant 
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uMd  for  isfiny  •xptrl— nf  1  iov^ttlgationt ,  ***  Stoker  [1937]  • 

It  attift  be  noted*  however*  thet  this  hydraulic  analogy  is  only 
valid  for  a  flat  horizontal  bottom,  as  may  be  seen  by  comparing 
(28.1)  and  (28.2),  and  even  wore  important  it  is  valid  only  as  a 
first  approziauition  even  for  the  nonlinear  case*  It  will  be  shown 
in  section  31  that  the  second  approxisatlon  to  shallow-w«ter  theory 
yields  finite  amplitude  waves  (the  sclicary  wave  or  cnoidal  waves) 
that  can  be  propagated  withcoit  a  change  in  shape  or  form,  a  fact 
which  completely  invalidates  the  hydraulic  analogy  to  compressible 
gas  flow  since  (28*2),  or  the  gas  dynamics  equation  ,  predicts  that 
any  finite  disturbance  quickly  forms  a  finite  discontinuity,  e.g. 

see  Lamb  [1932]  pages  278,  481 « 

In  section  29,  immediately  following,  it  will  be  shc^m  that 
even  for  the  linearized  first  approximation  the  hydraulic  analogy 
to  compressible  gas  flow  is  limited  to  a  flat  horizontal  bottom* 

29*  The  Linearized  Shallow-water  Theory* 

The  first  approximation  to  shallow-water  theory  c/m  now  be 
linearized  by  two  different  methods  that  are  suitable  for  various 
problems*  The  first  Mthod  is  more  appropriate  for  investigating 
steady  water  flow  in  canals  or  rivers  and  conelsts  of  the  following 
approximations  for  carrying  out  the  linearization 


u<x,«)  -  U  “  (29.1) 

'h  (x,*)  «  h(x,x)  (29.2) 

SO  that  (28,1)  is  linearized  to 

c 

c2  ••  gh(x;*)  (29.4) 


In  nrr - -  with  th*  previous  dlteustlon,  (29.3)  corrosponds  co  ch* 


linearized  gas  djnaMlca  equation  only  if  the  bottoe  is  flat  and 
horizontal  so  h  is  constant* 

The  second  Mthod  of  linearization  corresponds  to  the 
classical  t.  'al  wave  theory,  or  long  wave  theory,  e«g*  see  La»b 
[1932],  page  254,  and  can  be  obtained  by  writing 
u(x,*,t)  1.  I»(x,s,t)  <r<  1 

(x,z,t)  «  h(x,z) 

SO  that  (28.1)  is  linearized  to 


(29.5) 

(29.6) 

(29.7) 


(29.8) 


Again,  as  before,  (29*7)  corresponds  to  the  linearized  gas  dynamic 
case,  or  the  simple  acoustic  wave  propagation  equation,  only  if  the 
bottom  is  flat  and  horizontal.  In  this  case  the  general  solution 
of  (29.7)  for  one- dimensional  flow  is  the  wzll  known  d'Alembert 
solution  of  tha  simple  wave  equation. 

^(x,f;)  -  F(x-ct)  +  f(x4ct) 
c  -  rs  •  constant 
which  is  used  to  study  long  w^e  length  oscillations  in  canals 
when  the  water  is  either  at  rest  or  moving  with  a  velocity  of 
U  c.  The  limitation  to  linearized  perturbations  and  .ant  h 
for  one  dimensional  flow  allows  (28.1)  to  be  linearized  to 
^xt  -  -  -8  ^x 

1  «  "  -  J-ftt 


(29.9) 


and  vai'ioi&s  applications  of  this,  including  the  canal  theory  of 
tides,  are  given  in  Lamb  [1932]  pages  254-273. 

For  the  case  of  a  canal  having  a  non-rectanguiar  but  constant 
cross«sectiOD,  we  may  generalize  (29.9)  by  letting  h  be  the  mean 


dftpth  dftflxMd  4U  th«  undlstrubcd  cross -8«clIon«l  ares  (S)  dividad 
by  tha  width  (b)  of  tha  canal  at  tha  undi.stiirbad  fraa  water  surface, 
sea  Lai^  [1932]  page  256.  Whan  tha  canal  has  a  variable  depth  h(x) 
and  tha  disturbance  nay  be  considered  ona-diisansioital,  than  (29.7) 

■ay  be  written  In  tarss  of  the  varying  cross-sectional  area  S(x) 

1 

for  constant  width  b  as  follows 

S(x)  •  bh(x)i  b  «  constant 
Then  from  (29*9)  we  obtain 
^  x^x 

which  Is  tha  saaa  as  tha  expression  derived  by  Green  (1838)  for  a 
canal  that  Is  varying  In  both  width  b  and  depth  h  so  that 
S(x)  -  h(x)b(x) 

However  the  exact  linearised  first  order  approximation  is  (29.7), 
and  the  fora  of  this  equation  indicates  that  large  of  b'(x) 

would  invalidate  the  one -dimensional  ass\aiption,  especially  if 
is  relatively  large.  This  is  also  indicated  by  Land>  [1932]  page  274. 
However  (29.7)  provides  the  rigorous  proof  that  (29.10)  le  applicable 
to  one-dlisensiorial,  long  wave  lesigth,  small  amplitude  disturbances 
in  a  rectangulsr  cross-section  canal  having  a  constant  width  and  a 
varying  depth. 

Then  if  we  analyse  the  long  wave  lengths  having  a  simple 
harmonic  oscillation  with  a  frequency  of  a;/2T,  (29.11)  reduces  to 


tt 


(29.11) 
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If  we  aeeuiic,  folloirttig  Lamb  [1932]  page  275,  that  thla  expression 
gives  a  sablsfactory  approximation  to  a  varying  width  canal,  »*,an 
varies  linearly  with  x  ^lle  the  depth  h  remains  constant,  we 

2  I——  * 

obtan  Y  XX  +  i'yx  +(^)  y  -  0;  c  -  -  constant 


'uhose  general  solution  is  given  by  the  Bessel  functions  of  zero 
order  in  the  form 

^(x,t)  [AnCosttnt  +  BnSiii(4,t]  (29.13) 

^ere  in  most  applications  ®  since  -Yq(0)  -►  . 

On  the  other  hand  if  the  canal  width  remains  constant »  vdiile 
the  depth  varies  linearly  with  x  so  that  h(x)  -  h^x 
b  •  constant »  then  there  is  no  doubt  that  as  long  as  the  disturbance 
li  smell  and  one -dimensional  equation  (29.10)  is  applicable  and  may 
be  wrlttar.^  with  Cq  = /iSo  •  constant,  as 

which  has  a  solution  suitable  for  most  purposes  ss 

^(x,t)  ^  +  Dn^ol  lAflCoeWnt  +  Bn  sin  W^t]  (29.14) 


For  *noth*r  r.ctKigul.r  cross-scctlon  cast  let  h(x)^o(l-^), 

la* 

b  -  constant  so  that  (29.10),  with  -  1^  -  conatant,  bacomaa 
^  2 

V  o)/?  ^  ^  ■  ° 

la  «a 


2 

n(n+l)  •  (^L) 


Then  If  is  finite  for  /x/^  L,  a  solution  is  given  by  the 
Lsegendrs  polyncsdal  Pn(x/L),  of  degree  n  defining  the  natural 
oactlUtlon  £requ«ielat  by  -  n(n  +  1)  -  2,  6,  12,  20,  .  .  . 
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iVoaWnt  +  (29.15) 

Th«s«  last  problems  could  be  considered  es  the  long  period 
oscilletions  celled  seiches  which  occur  In  certain  lakes  or 
canals  throughout  the  world.  Many  applications  are  presented  by 
Chrystal  [1905»  1906]  and  the  periods  observed  in  several  lochs 
and  lakes  seem  to  correspond  to  those  calculated  by  the  linear 
shallow  water  theory.  The  linear  shallow  water  eqxiation  (29.7) 
should  be  very  suitable  for  the  study  of  seiches  because  of  their 
long  period  and  relatively  small  amplitude.  Usually  the  complete 
equation  (29.7)  must  be  solved  numerically  by  the  method  of  finite 
differences  because  the  contour  of  the  body  of  water  is  quite 
irregular  and  the  depth  variation  is  Important. 

There  are  several  simple  two-dimensional  cases  that  easily 
yield  explicit  solutions.  For  example »  the  long  period  oscilla¬ 
tions  in  a  rectangular  tank  having  a  constant  shallow  water  depth 
could  be  obtained  by  reducing  (29.7)  to 


where  the  prescribed  boundary  conditions  are 
^x(xo,z.t)  -  o  -  0  -^j(x,o,t) 

SO  the  appropriate  solution  for  ^  ox  ^  would  be  of  the  form 

2^  ^  Co.  (G^)  Co.  (n^  lA,^nCo.«Ji„nt)  +  J  (29.16) 

Similarly  for  a  circular  tank  having  a  constant  depth  (29.7)  may 
be  transformed  into  pooler  coordinates  (r^S  )  to  give 


The  solution  being  omitted  because  it  becomes  infinite  at  r  -  0. 
The  adnissible  values  of  (which  detemine  the  natural  frequen¬ 
cies  of  oscillation)  are  defined  by  the  boundary  condition 


0 


^  Ji 


or. 


for  n^O,  and 


0  or 


3.832,  7.016,  10.17, 


for  n  -  0  which  corresponds  to  the  syiametrical  oscillations. 

The  Yjj  solution,  %9hich  beliaves  as  (r^lnr)  as  r  approaches 
zero,  if  required  for  any  ring  type  body  of  shallow  water.  For 
example  an  island  in  the  center  of  a  circular  lake,  or  lagoon 
fomed  by  an  atoll  would  require  that  the  term  in  (29.17)  be 
replaced  by  the  complete  general  solution  of  Bessel's  equation  of 
ord.r  n,  nanely,  [Jn  +  EnYn  ^  n  rj  1  (29. IS) 


SO  as  to  be  able  to  satisfy  the  ring  type  boundary  conditions 


bj  +  bj  -  0 


Consequently  the  ring  type  boxmdary  conditions  are  satisfied  only 
by  the  particular  values  of  64},n  <^fined  by 

“^n  "  ^^n  (^(4B,n)  , 

TA(«  tii.n) 


(29.19) 
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Theo«  valuat  of  corrospond  to  the  neturel  frequencies  of  the 

long  period  oecilletion  of  the  ring  shaped  shallow  body  of  water* 

As  a  final  fxarsple  of  two-dinension^l  seiches  we  analyze  the 
long  period  siaple  harmonic  oscillation  in  a  shallow  circular  basin 
herring  a  depth  variation  dependent  only  upon  r  so  (29*7)  transforms 
into  polar  coordinates  as 

<hT  P  ^  ^  ^  °  (29.201 

Laab  [1932]  page  291|  has  shovm  that  the  solution  for  a  parabolic 
dapth  variation  h(r)  ■  hn  (1  -  £?.  )  i*  given  by 

^  (r.O.t)  -  F(eC.,4.r^>  (29.21) 

wh«r«  T  r*pr*s«nta  th*  hyp*rg*oa*trlc  series 


[  1  + 
a. 


(29.22) 


with  -i  -  n  (n-2)-m^ 

S^o 

c<o  1/2  (a  +  n),p  -  1/2  (2  +  m  -  n),^  -  (m  +  1) 

For  the  syaMtrical  modes  (m  -  0)  this  solution  rediices  to 
^  (r,t)  -  ^(Aj(  cos«„t  +  B],  ttnWnO  F(M,  l-N,  1>  ^  ) 

r  -  t  1-  «(«-!} (H-2}  Q  -  .  .  .  .) 

•2  .  4]!I(N-1)  -  8  ,  24  ,  48  ,  80  .  .  . 

8®o 

It  is  Interesting  to  compare  these  long  period  oscillations 
(^>>h),  wherein  the  Ilaid  particle  motion  is  approximately 
independent  of  the  depth  y,  with  the  short  period  oscillations 
(^<<  h)  which  have  a  direct  dependence  on  y*  For  example 
(23*14)  would  reduce  to  (29*16)  as  (y  +  h)— >0  since  coshy— >1 
and  tanhy  ,  »y  as  y-aO.  Also 

tsnK.  -VS  (l-^/^my  +  .  .  .  .  J  -yghll  +  o(iJ 


c 
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Similarly  (23.I5)  would  reduce  to  (29*  17)  and  (23*16)  would  be 
equivalent  to  the  additional  results  given  by  (29 .18)  and  (29*19)* 
This  eit^act  agreement  of  the  linearized  results  in  the  limiting 
case  is  encouraging  Justification  for  both  the  shallow-water 
approximation  and  the  infinitesimal -wave  approximation  since  they 
originate  not  only  from  different  physical  considerations »  but 
also  by  entirely  different  mathematical  Iterations  as  discussed 
in  section  10*  The  shallow-water  approximation  leads  to  hyperbolic 
type  nonlinear  equations »  while  the  infinitesimal -wave  approxima¬ 
tion  deals  with  linear  elliptic  equations.  Stoker  [1947]  page  32 
gives  a  detailed  comparison  of  the  two  linearized  approximations  for 
the  case  of  oaa-dimensional  flow  over  a  flat  bottom  at  a  6^  slope* 

29-  Linearized  Shallow-water  Tneory  Applied  to  Two-Dimensional 
Steady  Flow* 

The  first  method  of  linearizing  the  shallow-water  theory, 
as  given  by  (29*3),  would  be  applicable  to  the  determination  of 
the  variation  in  water  depth  for  the  steady  flow  in  a  shallow 
open  channel  or  river*  However,  in  practically  all  cases,  (29*3) 
must  be  solved  numerically  so  it  does  not  entail  a  prohibitive 
eniount  of  extra  labor  to  solve  the  more  exact  original  nonlinear 
flrftt-order  equation  (28*1)  directly  using  the  methods  discussed  in 
the  next  section  (30)  on  nonlinear  first-order  theory.  As  a 
matter  of  fact,  for  supercritical  flow  defined  by  U  >  |fgh,  the 
method  of  characteristics  is  very  easy  to  use  in  the  nuoMrical 
solution  of  the  flat  bottom  nearly  horizontal  open  channel  having 
a  varying  width,  as  shown  in  section  30*  The  subcritical  case, 
having  a  flow  velocity  everyi^ere  can  be  satisfactorily 
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«pproxlaAt«d  by  the  one- dimensional  hydraulic  theory  which 
aesuMS  that  the  velocity  at  each  cross-section  S(x)  is  indepen¬ 
dent  of  y  and  z*  This  method  of  course  would  yield  a  constant 
depth  over  a  given  cross-section  and  would  therefore  not  be 
satisfactory  for  predicting  the  rise  in  water  level  about  an 
island,  or  a  Jetty,  or  a  pile  in  a  swiftly  moving  relatively 
wide  stream*  For  this  particular  application  the  linearized  form 
of  (29.3)  is  very  useful,  especially  for  subcritical  flow. 

We  now  consider  the  application  of  (29.3)  to  the  problem  of 
determining  the  water  depth  variation  about  a  two-dimensional 
cylinder  that  is  perpendicular  to  the  bottom  and  has  a  narrow 


cross-section  parallel  to  the  flow  as  shown  in  Figure  38.  If 
the  bottom  is  approximately  flat  and  horizontal  everywhere  near 
the  vertical  cylinder,  then  we  may  consider  h  as  constant  and 
write  (29.3)  as 


/b  ^  4>jai  “  0  or  ^  -  0 

^  ^  “  (1*  -  constant  y  0 

The  ftmdamental  solution  of  (29.23),  in  view  of  (29.1),  for 
tifo-dimensional  profiles  that  may  be  considered  symoie^rical 
about  the  z-axis  as  shown  in  Figxtre  38,  is 

^(x,z)  “ li  +  (/S*)^ '  df 

V  r  ^ 


(29.23) 


X  - 


(X  -f  )2  +  (^z) 


(29.24) 


^  J  -}r  ^ 

Tl>e  boundary  condition  for  the  tvo^dlsensional  shape,  as  shown 
In  Figure  31,  is 
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(29.25) 


%9h«r€i  the  same  linearization  procedure  has  been  applied  to  the 
boundary  condition  as  was  used  in  deriving  (29.3).  'fherefore 
(29.24)  may  also  be  similarly  linearized  by  writing 

^  '4  /«|)  <^) 


sa  if  we  let  'x  •  ^ 


V/bZp  / 

p,  then  for  z^^  0 


f  -  A>  iM  ■'■  0(*^o) 

'  )  1  +  p"  2 

•A^.' 

Therefore  f(j)  -  (21J)  z^(^)  +  0(z2) 

SO  that  the  linearized  form  of  (29.24)  is 

'''*k  f - ?>  »;< 1. 1  + 


(29.26) 


x.z. 


jr/€> 


»  (x-/)^  +  (/2)Z) 


/  ^ 

w(x.z)  -  _z^  I  gp  (  ^  f 

^  ,  (x-  f)2  +(/S5 ) 


(29.27) 


^^z 

IT 


^  (x-f)^+^j)‘  IT 

On  the  actual  surface  of  the  two-dimensional  profile  (29.27)  may 
be  further  linearized  to  ^  ^ 

u(x,Zo)  ~  I  +  fhri  j  i-  f  1  . 


(29.28) 


w(x,z„)  -  zi(x) 


On  the  other  hand,  for  large  values  of  z  and  x«0  we  may  write 


(OtcP) 

a 

■  u 


T^3 


-C  , 


J/6  ^ 
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Th«  clumgc/h  (x,z)  In  the  original  constant  water  depth  h  can 
then  be  deteznined  by  the  linearized  relations  corresponding  to 


(29.1)  and  (29.2)  as 


'  <f> 

-■)  D 


(*.*) 


+  /-2l2g -  -  (^  +0(^)^(29.30) 

where  for  any  (x,z)  we  obtain and  from  (29.27).  The  term 
^  oust  be  included  in  (29.30),  althoiigh  must  be  omitted  in 
order  to  bti  consistent  with  linesrized  theory,  for  application  to 
large  values  of  z  since  (29.29)  shows  that  (^2^  whereas 

1/z^)  for  large  values  of  z.  However  for  small  values  of  z 
the  term  must  also  be  neglected.  For  example,  on  the  surface 


of  the  two-dimensional  profile  (z  -  Zq)  ,  (29.30)  reduces  to 

r  — I 


-  -  u2 


■  -t  +  "O-vJ  <«•”> 

siiice  and  are  both  of  0(Zq). 

These  relations  are  of  course  completely  restricted  to  flows 
that  are  everywhere  subcritical  since  (29.23)  shows  that  the 
Froude  nisiber  (F  -  U/  jT^)  must  be  everywhere  less  than  unity  to 
keep  0.  The  effect  of  increasing  Froude  number  is  to  increase 
and  therefore  decrease^  ,  since ^  decreases.  It  is  seen  that 
this  effect  increases  as  z  increases,  the  greatest  effect  being  on 
^  ^  limiting  case  of  very  large  values  of  z  as 

shown  in  (29.29).  This  relation,  or  preferably  (29.27),  could  be 


used  to  predict  the  additional  change  in  /^(ic,y)  d\ie  to  a  finite 

stream  width  by  using  the  incitement  of  from  one  mirror  image  to 

represent  die  first  approximation  to  the  channel  boundary  wall  as 

indicated  in  Figure  38.  For  slender  cylinders  in  a  narrow  channel 
the  "one-dimensional"  approximation  of  section  30 -/^is  generally  used, 
this  allows  an  approximation  for  frictional  head  loss  which  becomes 

relatively  more  important  as  the  channel  width  decreases. 


k 


E-14-A 


For  supercritical  flow  (F  «  U/( gh  (29.23)  must  be 


ir^ltCen  as 

-<l>zz  ' 
-  (F^  -  1) 


or 


£l- 

'd  / 

1)  -  constant  0 


(29.32) 


Wow,  however  (29.32)  cannot  provide  a  satisfactory  approximation 
of  the  change  in  water  depth  at  some  distance  from  the  two-dimensional 
profile  since  its  general  solution  is 

^(x,z)  -  G(x-Bz)  g(x  +  Bz)  (29.33) 

which  predicts  no  change,  even  upon  approaching  infinity,  along  the 
lines  of  constant  slope  dz/dx  -  ±  ^  ■  ± 

Consequently  the  nonlinear  method  of  characteristics,  as  will  be 
described  in  section  30,  must  be  used  in  predicting  the  depth 
'Viiriation  at  any  finite  distance  from  the  profile.  Although  the 
method  of  characteristics  will  directly  and  easily  give  the  velocity 
distribution  or  depth  variation  on  the  profile  itself,  we  will  now 
derive  the  linearized  variation  valid  on  the  profile  surface.  The 
result  will  be  of  crucial  importance  in  evaluating  the  validity  of 
the  nonlinear  first-order  shallow-water  theory  (28.1),  since  any 
discrepancy  between  the  linearized  result  and  the  nonlinear 
results  from  (28.1)  would  indicate  that  the  perturbations  Involved 
are  sufficiently  large  that  the  second-order  shallow  water  theory  of 
section  31  must  be  introduced. 

The  linearized  solution  of  (29.32)  for  any  sharp-nose  slender 

two -dimensional  profile,  as  in  Figure  36,  is  obtained  from  the 

general  solution  (29.33)  and  the  linearized  boundary  condition 

(29.24)  as  follows:  Zq(x)  -  ^z(x,Zo)  -  -BG*(x-Bzo);  z-Xq  ^  ^ 

ij  IT 
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■  I 

thsrafor*  G  (x-Bzq)  "  “  y 

MO  that  on  the  profile  surface  (z  «  Zq) 

u(x,Zo)  -  U  +  +  G*(x-Bzo)  -  Ull  -  2q(x)  ] 

'  “B 


w(x,Zq)  - 

Then  the  variation  In  water 
by  (29.30)  at 


(29.34) 

depth  on  the  profile  surface  la  given 


^  [  *o<*)  +  0(Zo)^J  (29.35) 

For  flow  that  is  everywhere  supercritical  so  B^  -  (F^»  1)?0 

It  should  be  noted  that  (29*23)  and  (29.32)  are  identical  to 
the  linearized  potential  equations  for  two-dimensional  steady  sub¬ 
sonic  flow  and  supersonic  flow  respectively  if  we  simply  replace 
the  Froude  number  (F  -  U/j[^)  by  the  Mach  ninber  (M  -  U/C), 
see  (28.3).  This  is  in  complete  accord  with  the  statezaent  that  the 
hydraulic  analogy  is  valid  for  the  flow  over  a  flat  horizontal  bottom 
(which  is  equivalent  to  the  two-dimensional  isentropic  flow  of  a 
fictltous  perfect  gas  having  a  specific  heat  ratio  of/'  -  2). 
Consequently  equations  (29.24)  through  (29.29)  are  Identical  to  the 
subsonic  flow  about  slender  two-dimensional  profiles  in  free  air  or 
In  a  rectangular  cross-section  wind  tunnel  as  derived  by  Laitone 
( 1946] .  These  equations  confi  rm  the  known  result  that  the  linearized 
•quations  are  independent  of  the  value  of  the  specific  heat  ratio  , 
Similarly  equation  (29.34)  Is  identical  to  the  well-known  linearized 
tnro-dimansional  supersonic  flow  solution  if  we  let 
-  1)  ■  (m2  -  1) >  0. 
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Althoxigh  these  linearized  results  are  very  satisfactory  for 
slender  sharp  nose  profiles,  they  only  apply  for  Froude  numbers 
thMt  are  not  too  near  unity.  That  is  they  are  not  applicable  to 
flows  near  critical  (U  c,  equivalent  to  sonic  flow). 

For  these  cases  we  must  return  to  the  nonlinear  equation  (28.1) 
as  discussed  in  section  30. 
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30  “  Nonlinear  Shallcw-water  Theory. 

Thl^  section  will  primarily  discuss  methods  for  obtaining 
solutions  of  the  nonlinear  equation  (28.1)  which  provides  the  flrst- 
order  approximation  of  the  shaliow-water  theory.  The  special  cases 
to  be  considered  are  the  one-diaaenslonal  unsteady  flow  and  the  two- 
dimarslonal  steady  flow  In  open  channels.  This  will  provide  a  basis 
for  dliicusslng  the  one- dimensional  assumption  of  open  channel  flow. 
Finally  the  hydraulic  Jumps,  and  their  relation  to  the  first-order 
shallow-watar  theory,  will  be  discussed. 

yO  ~  One-dlmenslonai  Non-steady,  First-order,  Shallow-water  Theory. 

By  Agfualng  one-dlmensloqal  flow  In  the  x  direction  only,  the 
nonlinear  equation  (28.1)  redxices  to 

Ut  +  uu^  +  g  (^  +  -  gh^  (30.1) 

(/I  +  h)t  +  M*]+  h)l^  -  ht  -  0 

again  It  should  be  noted  that  these  are  equivalent  to  the  gas 
dynamics  equations  ujwjn  Introducing  (28.3),  only  If  Che  bottom  Is 
flat  and  horlsontal  so  h^  “  0. 

Now  If  we  let 

c2(x,t)  -  g  (  Y  (x,t)  +  h(x)  ] 

2cCjj  -  8(^7+ h)jj  (30.2) 

2ccj;  ■  +  h)^ 

and  give  the  Initial  conditions  as  du/do<^  and  dc/d^x  along  a  space 
curve  defined  by  x<o<),  t(o^),  then  we  may  write  (30.1)  as 
u(ux)  +  (Ut)  +  2c(c,)  +  0  -  gh^  ;  c(u^)  +  0  +  2u(cx)  +  2(c^)  -  0 

<«x>  +0  +  ° 

0  +  0  +  x^  (c,)  +  t,  (Ct)  - 


(30.3) 
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This  set  of  four  equations  can  be  solved  uniquely  for  u^,  u^., 

^x>  long  as  the  determinant  of  the  coefficients  In  (30.3) 

does  not  vanish.  Consequently  this  necessary  condition  Is  violated 
along  the  characteristic  curves  x(o<),  t(o<)  defined  by 


(30.4) 


vblch  may  be  easily  expanded  by  the  minors  of  the  bbttocn  row  to  obtain 

+  (u2  -c2)tp,  2  ]  . 

-  -  (u  -  c)t^  1  [x^  -  (u  +  c)t^  ]  -  0 

Therefore  the  characteristic  curves  are  defined  by 

“  (^)  “  l'*^**^^  1  c(x.t)  1  (30.5) 

Since  hx  Is  given,  and  appears  only  on  the  right  hand  side  of  the 


first  equation  In  (30.3)  therefore  the  characteristic  curves  as 

i 

defined  In  (30.5)  are  Identical  to  those  In  the  gas  dynamics  case, 


e.g.  see  Courant  and  Friedrichs  [1948].  However,  the  Rlemann 

t 

Invariants,  or  quantities  that  can  be  constant  along  a  characteristic 
curve,  now  depend  upon  the  bottom  slope^  as  may  be  seen  by  adding  the 
tvo  equations  In  (30.1)  after  Introducing  (30.2)  so  as  to  obtain 


(u  +  2c)^  +  (u  +  c)  (u  +  2c)^  -  gh^  - 
"  ^  ^  +  (u  +  c)-^  ]  [u  (x,t)  +  2c(x,t)  ] 


(30.6) 


These  give  the  same  Rlemann  Invariants  as  in  the  Isentroplc  one- 
dlaienslonal  unsteady  gas  flow  with  a  specific  heat  ratio  of  2 


only  If  h^  "  0,  e.g.  see  Courant  and  Friedrichs  [1948]  page  87.  No 
elaple  Rlemann  Invariant  Involving  only  u  and  c  Is  possible  If  h^ 
varies  with  x,  however  If  h^  Is  constant,  so  gh^^  ■  n  ■  constant. 
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th«n  (30.6)  m£y  b%  written 

'  7%  +  (u  +  c)  -^  ]  [  u  +  2c  -  mt)  -  0  (30.7) 

Slailerly,  by  subtracting  the  two  eqiietions  in  (30.1),  we  obtain 

[  -T%  +  (u  -  c)  ,2^  ]  [u  -  2c  -  »t]  «  0  (30.8) 

Consequently  the  basic  statements  relating  tho  characteristic 
curves  and  Riemann  invariants  of  equation  (30.1)  with  gh^  ^  m  «  constant 
oay  be  auoaiarised  as  follows: 

(u  +  2c  -  mt)  ■  R(x,  t)  •  constant  along  a  curve 

defined  by  -  (u  -f- c)  ^3^  ,^ 


(u  -  2c  “  mt)  -  -SCxjt)  constant  along  a  curve  C_ 

defined  by  ^  "  (u  -  c) 

Figure  39  shows  typical  sets  of  curves  in  the  (x,t)  plane.  The  above 
equations  show  that  in  any  given  region  in  the  (x,t)  plane  there  are 
three  basic  types  of  solutions,  namely: 

(1)  The  constant  steady  state  in  which  u  and  c  remain  constant 
averyuhere  in  the  region  so  all  characteristics  form  straight  lines, 

(2)  The  general  flow  in  %dtich  neither  R  nor  S  is  constant  in 
a  finite  region, 

(3)  The  special  c^ee  of  a  simple  wave  over  a  flat  horizontal  bottom 
(m  "  0)  wherein  a  constant  steady  state  region  is  separated 
from  a  varying  region  by  a  straiglit  characteristic  line  along 
which  either  R  or  S,  is  constant. 

The  first  type  solution  obviously  has  R  and  S  constant  throughout  the 
region  only  if  the  bottom  is  flat  and  horizontal  (m  «  0)  •  The  second 
type  of  solution  is  coRplicated  and  can  best  be  obtained  by  the  method 
of  finite  differences,  e.  g.  see  Stoker  [1957]  pages  293-300.  The 
third  type  of  solution  will  now  be  discussed  since  it  has  considerable 
physical  significaz>ce  for  many  probloms  concerning  the  propagation  of 


i 
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A  disturbATice  into  water  that  is  originally  at  constant  depth  and 
con#t:Ant  velocity,  and  extends  an  unlimited  distance  for  x  >  0. 

When  a  disturbance  moves  into  still  water  at  constant  depth 

over  s  flat  horizontal  bottom  (m  «  0) ,  then  it  is  obvious  that 

((ix/dt)^  -  c(«o)  is  the  characteristic,  which  is  now  a  straight 

line  9  that  vmst  continually  separate  the  steady  state  region  from 

the  disturbance  region  in  the  (x,t)  plane  as  indicated  in  Figure 

39.  This  characteristic  curve  must  be  a  straight  line  since  there 

Is  a  constant  steady  state  always  ahead  of  it  so  that  (dx/dt)^  ■ 

constant  and  therefore  x^  ■  c(c>^)t.  Also  either  R  or  S  must  be 

constant  along  the  characteristic,  and  since  corresponds  to 

or  (dx/dt)^  «c(®^)  >0,  as  in  Figure  39,  therefore  R^  -  2c(cio) 

■  constant.  This  type  of  simple  wave  having  (—1  ■  c(o^)>  0 

'dt/o 

and  ••  2c(cio)  -  constant  is  called  a  forward- facing  wave  since  the 
partlcL^  paths  enter  from  the  side  with  greater  values  of  x,  as  in 
Figure  39.  The  value  of  R  varies  on  each  C4.  characteristic  inside 
the  regicn  of  the  disturbance  since  u  and  c  both  vary  due  to  the 
disturbance  and  none  of  the  characteristic  lines  can  evex*  inter¬ 
sect  CJ.  However,  every  characteristic  terminates  at  CJ,  as  shown 
in  Figure  39,  and  since  S  remains  constant  on  any  given  C.  character¬ 
istic  curve,  therefore  S  is  everywhere  constant  since  every  C. 
characteristic  must  have  the  same  value  S(x,t)  R^  ■■  2c{oCj  m  constant 

/tO 

on  , 

The  same  considerations  are  true  even  if  the  steady  state  constant 
depth  water  into  which  the  disturbance  is  being  prop^ated  is  flowing 
with  a  constant  velocity  u(^)  The  only  change  is  that  now 

the  following  are  constant, 


E-21 


f^]  -  (u(oo)  +  c(oo)]  >0,  R„  -  l2c(o<=)  +  u(oo))  on  Q®  only, 

^dt/o 

while  on  ell  C.,  [2c(x,t)  -  u(x,t)]  -  S(x,t)  -  [2c(o^’)  -  u(oo)l  - 
constsnt.  Siailarly  all  R  In  the  disturbance  region  vary  as 
R(x,t)  -  (2c(x,t)  +  u(x,t)] 

as  indicated  In  Figure  39  for  the  simple  forward- facing  (C^  wave. 

Xs  shown  in  Flgvire  40. a  simple  backward- faclng(C°) occurs  if  R-coiistant  _ 
T1ieS€  WAV6S  Ar®  called  simple  waves  because  all  the  chaicactei^j. sties 
of  the  family  which  has  the  Riematm  invariant  take  on  a  different 
cunitant  for  each  line^  must  form  straight  lines.  For  example, 
referring  to  Figure  39,  the  forward- facing  waves  (dx/dt>0)  have 
S(x,y)  constant  everywhere  and  R(x,y)  varying  so  the  characteristics 
form  straight  lines.  On  the  other  hand  in  Figure  40,  the  backward- 
facing  wave  (dx/dt<r0)  has  R(x,y)  constant  everywhere  and  S(x,y) 
varying  so  now  only  the  C  characteristics  form  straight  lines.  The 

teristics  of  one  family  only  must  form  straight  lines  in  a  simple 
wave  because  only  one  of  the  Riemann  invariants  (S  or  R)  is  constant 
in  the  entire  region  of  the  disturbance.  For  example,  in  the  case  of 
the  forwaid- facing  simple  wave  in  Figure  39,  we  have  S  constant  in 
the  region  of  the  disturbance.  Therefore  from  (30.9)  and  Figure 


39  we  may  write, 

“Si  "  “S2  “  ''^'3  “  “S4  -  (ui  -  2ci)-(u2-2c2)"(u3-2c2)-(u4-2c4)-constant. 

Rj^*R3“(ui  +  2ci)«(u3+2c3)  R2  “  ^4  “  (^*2  ^02)  ■  (^4  +  204) 

Consequently,  u^  -  U3,  c:  •  C3,  U2  “  C2  "  C4,  and 


ui  ^  U2,  Cl  ^  C2>  ^  U4,  C3  ^  C4  so  that 


(ui  -  C|)  r  (u2  -  C2)  ^  ccniitant,  so  o.  is  curved 
m  -f  c^)  *  (u3  +  C3)  "  constant,  so  is  straight 


(U2  +  C2)  -  (u4  -(•  C4)  -  constant  so  is  straight. 
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It  is  Important  to  note  that  these  simple  waves  can  exist  only 
over  a  flat  horizontal  bottom  so  m  -  0. 

We  have  now  shown  how  the  method  of  characteristics  for  one- 
dimenslonal  unsteady  flow  has  resulted  in  the  concept  of  the  simple 
wave  which  quickly  gives  a  numerical  evaluation  of  the  propagation 
of  a  one -dimensional  disturbance  into  water  at  a  constant  depth 
that  is  moving  at  constant  speed.  The  solution  of  this  problem 
in  the  (x,t)  plane  can  be  obtained  by  direct:  application  of  (30.9). 

For  example,  the  usual  case  of  a  forward- facing  wave  having  S  every¬ 
where  constant,  and  straight  C4.  characteristic  lines,  as  shown  in 
Figure  39,  has  the  slope  of  the  C^.  straight  lines  determined  directly 
by  the  time  history  of  the  disturbance  at  x  ■  0,  and  equations 

(30.2)  and  (30.9)  which  show  that 

»  u(0,t)  +  c(0,t)  •  constant  -  u(0,t)  +yg[h  +  (0,t)]  (30.10) 

and  along  any  given  straight  line  having  this  constant  slope 

R(3c,t)Q  -  u(0,t)  +  2c(0,t)  -  constant  -  u(0,t)  +  1  V  g[h+  ^(0,t)] 

^  ...  (30.11) 
Consequently  the  values  of  u  and  c  are  deteiuilned  in  the 

entire  region  shown  in  Figure  39  by  the  given  values  on  the  t-axls. 
The  curved  C_  characteristics  need  not  be  calculated,  since  the 
desired  niznerical  solution  is  independent  of  them.  Their  existence 
however  can  lead  to  a  simplification  in  the  nixnerical  calculation  of 
(30.10)  since,  in  the  case  shown  in  Figure  39,  each  curved  char¬ 
acteristic  extends  from  the  C°  characteristic  to  the  t-axis,  and  on 
each  and  every  C.,  characteristic,  -S  «  [u(cxo)  -  2c(oo)]  ■  constant. 
Therefore  at  every  point  on  the  t-axls  that  can  be  reached  by  a  C_ 
characteristic  we  must  have 

S  -  [2c(0.t)-u(0,t)l  -  [2c(oo)-u(oo)]  -  [2  -  u(oo)]  -  constant 

(30.12) 


E-23 


Of  coursa  the  c_  characteristics  cen  continue  from  to  the 
t-axle  only  if  or  u<c,  so  that  in  this 

case  (30.10)  may  be  simplified  to 
/^]  ■  u(0,t)  +  c(0|t)  -  constant  - 

(30.13) 

-/C’3/2)u(0,t)  -(i/^tu(o5)-2c(oo)]J  -  f(3/^(0,t)Hi/^(oo)4-l/lh  |- 
-|  3c(0,t)  +  [uCoo)  -  2c(oo)]j  ••  [3\/gth  +'^<0,t)l '  +u(oo)-2  ? 

Consequently  the  problem  is  solved  in  the  region  so  defined  if 

I 

either  u(0,t)  or  c(0,t)  is  alone  given.  The  surface  elevation 


would  be  given  by  (30.2)  as 

h+7(x,t)  .  ,  h  -  -  conatTOt 

‘  g  '  O 


(30.14) 


in  every  case  of  disturbance  propogationa  into  a  constant  water  depth 


over  a  flat  horizontal  bottom  so  m  -■  0. 


Many  other  physical  problems  can  be  simulated  by  giving  the  data 
along  a  prescribed  curve  (x,t)  for  x  -  0,  e.g.,  see  Stoker  [1957] 
where  the  disturbance  created  by  the  breaking  of  a  dam»  and  the  effect 
of  moving  a  vertical  end  plate  in  a  rectan^^ular  cross-section  tank 
of  still  water,  u(''^)  -  0,  are  considered.  Since  the  bottom  is  flat 
and  horizontal  and  m  -  0,  all  of  the  equations  following  (30.9)  are 
equivalent  to  the  gas  dynamics  equations  xfith  a  specific  heat  ratio 
of  2,  conaeqviently  the  problems  in  finite  channel  lengths  which 
produce  wave  reflections  at  either  end;  as  solved  in  Courant  and 
Friedrichs  [1943],  are  also  applicable.  In  this  hydraulic  analogy 
to  compressible  flow  it  is  important  to  note  that  (30.13)  is  only 
applicable  to  subcritical  flow,  which  is  equivalent  to  subsonic  gas 
flow,  since  we  must  have  "  (u-c)<^0,  or  u(oo)<c(oo)  -  ^ST. 


i 
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When  the  flow  Is  supercritical,  so  that  u(oo)  >  c(oo) 
corrsspcndlng  to  supersonic  gas  flow,  then  the  slope  of  both  the 
4ind  C„  characteristics  are  negative.  Consequently  the  two 
families  can  meet  in  a  cusp,  and  the  C„  characteristics  cannot 
intersect  both  the  t-axis  and  the  undisturbed  steady  supercritical 

I 

state  that  lies  at,  and  to  the  right  of,  C^.  Therefore,  in  order 

1  V  ^  ^0  ..  1  ^  ^  Cl  •  ^  ^  «  1  g%  ^  ^  ^ 

a.xwiiv  ^&a5  iJL,  UliCS  CwilOL^lU 


value  of  S,  as  given  by  (30,12),  must  be  very  carefully  defined. 

Another  limitation  on  all  the  preceding  equations  is  indicated 
for  the  compression  wave  depicted  in  Figures  39  and  40.  This 

t  *  ^ 

limitation  is  defined  by  the  envelope  of  the  straight  characteristic 
lines  that  must  always  form  for  a  compression  wave  in  this  first 
order  theory,  as  will  be  proven  later.  This  envelope  of  the  straight 
characteristic  lines  corresponds  to  a  discontinuity  that  can  be 
interpreted  as  a  discontinuity  in^,  or  the  breaking  of  the  wave 
crest.  This  leads  to  the  hydraulic  Jump  or  surge  that  will  be  discussed 
lat;er.  The  gas  dynamic  case  has  the  envelope  of  the  straight  character¬ 
istic  lines  interpreted  as  a  steady-state  shock  wave,  e.g.  see 
Courant  and  Friedrichs  [1948]  pages  110-181. 

30-^  -  Two-dimensional,  Steady,  Supercritical  FIca^  by  the 
First-order  Shallow-water  Theory. 


We  will  now  investigate  the  characteristic  curves  of  the  nonii.iear 
tquations  of  the  first-order  shallow-water  theory  for  the  case  of 
sceady  two-dimensional  flow.  We  will  find  that  real  characteristic 
curves,  that  are  a  great  aid  to  numerical  calculations,  exist  only 
In  the  regions  wherein  the  flow  is  everywhere  supercritical. 


If  we  consider  the  steady  two-dimensional  flow  over  a  flat 
horizontal  bottom,  then  we  may  write  (28 •2.)  as 
uu^  +  wu^  -  +  ho)x  "  -(c^)x 

uw^  +  wWj  “  (30.15) 

+  ho)]^  +  [w('^+  ho)]^  -  0  ■  (uc2)^  +  (wc^)^,  “  0 

u  -  0^,  w  -  0j,  Uj  -  Wx  “  0XZ 


By  multiplying  the  first  equation  by  u  -  0^  second  by  w  -  0^, 

and  adding,  we  obtain 

[0^0^  +  20^0x0„+  0^0„  1  -  -  I0^(c2)^+  9l(ch^]  - 

-  ( C^xx  +  ^2z) 

Therefore 


+  2 


(30.17) 


or, 

(1-^  )0xx  -  2  ^  0X2  +  (1  -  f?  ><»**  “  0 

where  c^(x,z)  -  +  ^(x,z)]  and  h^  now  is  the  still  water  aepth 

f  jwd  whenever  (u^  +  w^)  -  0  -  ^  .  Note  that  (30.18)  immediately 
linearizes  to  (29.3),  so  the  numerical  differences  between  the 
solutions  of  (29.3)  and  (30.18)  will  provide  an  estimate  of  whether 
or  not  the  second-order  shallow-water  theory,  as  discussed  in  section 
31,  must  be  introduced. 

The  characteristic  curves  of  (30.18)  may  be  found  in  a  manner 
similar  to  that  used  for  (30.3)  by  finding  the  space  curve  [x(<?<), 
^(c^))  along  which  prescribed  values  of  0^^  and  0^  cannot  determine 

®xx*  ^xz*  ^zz* 


Therefore  we  write 
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)^xx  )^xz  ^ 

a  <x 

0  +  (Xo<)^xz  +(^of)»zz* 

dor 

which  may  not  hcive  a  solution  if  the 
zero,  that  is  if 

(1  .  ^  )  (  -  2uw  )  (  1  - 


(30.19) 


determinant  of  the  coefficient  Is 


(30.20) 


(30.21) 


which  therefore  givas  the  slope  of  the  two  families  (c^  and  c_)  of 

characteristic  curv'es.  Now  however,  entirely  unlike  the  previous 

one -dimensional  unsteady  flow  solution,  the  characteristics  curv€:s 

2  2  2 

exist  only  for  supercritical  flow  so  that  (u  +  w  )  ^  c  ■  g(hQ+  J;  ) . 
The  fact  that  the  characteristic  curves  are  real  for  supercritical 
flow  means  that  in  this  case  the  nonlinear  eqxiatlon  (30.17)  is 
hyperbolic.  However  for  subcrltical  flow,  since  the  characteristic 
curves  are  then  Imaginary  (complex)  functions,  it  is  of  the  elliptic 
type,  e.g.  see  Courant  and  Friedrichs  [1948]  pages  40-55,  or 
Preiswerk  [ 1938] . 

We  can  obtain  a  solution  for  the  variation  of  the  quantity 
rix.z)  -  ^  1  (30.22) 


V  w 


a 


•A/ 

(Which  cWtfinaa  the  Prouda  nvnber  of  the  tiipercrltical  flow)  along 
charts terii tic  curve  by  tranifonaing  (30.18)  into  tha  hodograph  (u,v) 
plana  through  tha  uaa  of  tha  Lagandra  coiic«ct  transformation  which  is 
given  by,e.g.  tfe  Courant  and  Priadrichi  [19481  PAg«  249  or  Praiswark 


[1938]^  (x0x  +  *02  "  0)  ■  (xu  +  2w  -  0) 

d  ^  •  (xdu  +  udx  +  xdw  +  wdz  -  d0)  •  (xdu  +  zdw) 
th.refor.i  X  -  ^  y,  *  dx  -  j^du  +  \ 

dy  -  z^da  +  ^di»  u»,du  +  ^ 

Solving  for  du  and  dw  wa  obtain 

du  -  l/M  (^wwdx  -  dS»3j  -  (035xdx  +  9^dz) 

dw  -  1/M  (  -^uwdx  +  ^uud*)  -  d0z  -  (^xz**  +  ^zz***) 


whara  H  * 


^  uu  "3f 


uw 


•o  that 


^uw  1(ww 

a  • 

^XX  11 


*  0 


^xs 


-  .^kipc* 


zz 


and  tharafora  tha  nonlinear  aquation  (30.17)  in  tha  physical  (x»s) 
plane  is  transfomad  into  a  linear  aquation  in  tha  hodograph  (u,w) 
plana,  as  givan  by 

(£.  -  l)\  -au  -  2  3  ^  uw  +  °  <3®-23) 

c  ^ 

Th.  a«ie  proc.dura  ms  ua.d  in  (30.19)  through  (30.21)  or  «  almple 
cooparUon  of  (30.17),  (30.21),  «td  (30.23),  .how.  that  the  character¬ 
istic  curves  of  (30.23)  In  th.  hodograph  (u,w)  plan,  are  defined  by 


lii 

u 


-S'- 


J 


(30.24) 


Tha  charactariatic  curva  P-  in  tha  hodograph  (u,t/)  plane  is 
orthogonal  to  tha  charactarintic  curva  C+  in  the  physical  (x,2) 
plana  if  wa  suparlmposa  tha  two  planes  so  that  tha  velocity  vectors 
coincide.  This  may  be  aaaily  shown  by  rotating  tha  axes  for  (30,2 l)dr>dl 
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(30.24)  80  that  w  -  0  (sae  Figure  41  ),  end  the  equations  for  the 
slopes  of  the  cherecterlstlc  curves  end  p-  simplify  to 

S'  -  ‘ 


c+ 


p2  -  1 


•  ) 

du  'n_ 


(30.25) 


Similarly  /  +  is  orthogonal  to  C_  ^en  the  planes  are  superimposed 
so  the  velocity  vectors  arc  coincident  (see  Figure  41). 

Equations  (30.24)  and  (30.25)  show  that  along  any  characteristic 
curvs  there  exists  a  simple  solution  which  Is  Independent  of  tne 

boundary  conditions  of  a  particular  problem,  since  we  can  directly 
Integrate  (30.25) 


for  axes  rotated  so  w  -  0  (dw  -  ud®),  see  Figure  41.  We  thereby 
integrate*  (30.25)  as 


dS 


r- 


r- 

\A  «!■  [/T 


3“ 


Con8€qu*ntly  (30.26)  provides  a 


tan'^  Vf^-i'j  sr  f(^P)  (30.26) 

solution,  independent  of 


the  boundary  conditions  in  the  physical  plane,  for  any  two-dimen¬ 
sional  potentiil  flai#  that  possesses  the  property  of  having  simple 


wavea  in  the  given  region  so  that  the  velocity  vector  follows/" 


In  the  hodograph  plane.  The  nxnerical  values  from  (30.26) 


*  see  (30.27)  ai'.d  (30.29)  which  show  that  with  w  -  0 


are 
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1.11  Flfurtt  41  Axid  iirfi  t^bulfttsd  in  Tftblft  I  (t®k®n  ri'oiii 
Fr«lMtrk  [1938]. 

Th«  ttstful  r«litlon  b«twt«n  c  «nd  (u,w)  that  w*8  used  to 
Inttgrtt*  (30.26)  and  calculate  Table  I  Is  obtained  by  ultlplying 
th«  flret  equation  In  (30  *15)  by  dx  and  the  second  ec  ^tlon  by  dz, 
and  then  adding  then  so  as  to  obtain 

uCu^jdjj  +  *♦*  <c^)2dz] 

udu  4-  wde  •  -  d(e^)  *  -  gd>^ 

1/2  (u^  4  u^)  +  c^  -  constant  ■  l/2<u^  +  w^)  +  s(^o  ^  )  (30.27) 


Therefore 

(nL+jd)  +  c2  -  +  g<ho  +'yj  )  - 


1^-2  I 

)■=  =  df. 


a 

9 

max 


<30.28) 


■  <8ho)  “  “  <1  ® 

«ih«r.  <»•«  ngur.  42)  ho  t«  the  still  weter  depth  (or  stagnation 

total  head  depth)  that  corresponds  to  (u^  +  w“)  -  0  - 
(u2  +  t»^),  is  the  llaiting  resultant  velocity  squared  that  is 
approached  .hen  the  flowing  water  depth  approaches  zero  -ho), 

and  is  the  critical  <F-1)  speed  when  the  resultant  velocity 


4  W 


eb  f 


is  critical  <F-1),  so  that 


"  (  8^**o 


-  ]f?2/^gho)'  -  ^2/^  +  g(h^+ 


(u^  +  w^) 


/ 


_  ®-K-  -  2/3  -  -2 

gho 


(!»/ .  -  (1 .  ^ -  (1  *  5^)  -  <|»,)'<^) 


3 


<3/^ 


(30.29) 

ho 


p2  . 


3  -  V 

f  3r^ 

I  2+F^ 

•  «« 


(30.29) 
cont ' d 


It  IB  very  ut.ful  to  not*  that  equations  (30.26)  through  (30.25)  mey 
*11  be  obtained  directly  froo  the  two-dimenalonal  laentroplc  gas 
flow  equation*  by  slaply  letting  the  specific  heat  ratio  -  2 
«d  r  =  M,  E  =  ;  e.g.  compare  with  Cour«»t  and  Friedrich*  119481. 

as  had  previously  shown  by  Prelswerk  [1938]. 

The  Rlemaim  Invariants  for  the  characteristic  curves  (C+,  C.) 
will  now  be  determined.  First  we  can  show  that  the  velocity  component 
normal  to  the  characteristic  curve*  Is  always  the  local  velocity  of 
the  shallow-water  wave  propagation,  c(x,z),  by  writing  (30.21)  a* 

(udz  -  ■  e2j(dx)2  +  (d*)^3(j^-  c^(dA)(;^ 

since  the  relation  between  the  normal  direction  (n) ,  and  the  tangen¬ 
tial  direction  (A)  along  the  characteristic  curve  (C+)  U  given 

b,  (».  ngur.  41)  m  -I'frL  .  1  - 


b,  (».  ngur.  41)  -  -<W)c. 

Similarly,  If/l  1*  the  tangential  direction  along  C., 
-  (u^  ■ 

Also,  from  Figure  41  and  (30.21) 

tan  -  (fes.o  -  +  ,  •!««  -  ■  q 


■  0 

Ci^.^ 


f  I  ■  ^  £ 


(30.31) 


where  q  Is  the  resultant  velocity  magnitude  so 
^2  .  .  (j^  +  c2)  -  (9^  +  c2)j  e  - 

u  -  q  CCS  9,  w  -  q  sin  0, 

9^  -  q  cos  0<  -  ^ 


(30.32) 
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Subitttuting  (30.31)  and  (30.32)  Into  (30.21)  and  (30.24)  wa  obtain 


Th«r«forft»  proven  b«for«  in  (30.25) » 


•  tan(0  !■  )  (30.33) 


■  -  cot(0  ±  ex' )  (30.34) 


(30.35) 


that  i»,  a.  .how»  In  Fijuri  41,  the  C+  charact.rlstlc  curv.a  In  the 
phyilcal  ^x,x)  plana  ar.  at  avary  corraaponding  point  orthogonal  to 
the  r.  oharactarlatic  curv.a  In  the  hodograph  <u,w)  plan*.  All 
th.aa  raaulta  are  the  aaaa  aa  In  the  gaa  dynamlca  caaa  »ihar.  tha  C+, 
charactarlatlc  curv.a  are  referred  to  aa  the  Mach  linea  since,  aa 
shown  by  (30.30),  the  normal  velocity  component  la  always  the  local 


•p««d  of  lounde 

Now,  as  shown  in  Figure  41, 


(^)p  ■  ^ 


«n  oOp  ;  U)^  -  (0  +  ex'  -  -J-) 


tan  uj-  ; 


(Q  -  (X  +  --J-) 


therefore  (30.33)  may  be  written  as 


..  (9  +  o<' )  -  tan  ( tOp  +  ) 

{ “  tan  (0  -  o<  )  “  tan  ( ) 

Consequently  tha  Ilemann  invariants  are  given  by 
R  -  (9 +o<*  60^  • -5^  )  S  -  (0  -CK  -  COp^ ) 
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This  ns7  b« 


from  (30.24)  and  Figure  41  since 

j  cano<  I  ■  f+l  cot 


or  scs  Coursnt  and  Friedrichs  (1948]  page  266.  - - - 

Thsre£or«  g  •  ^  -I-  tan"^  ****  '  3*  ' 


(IpT .  ,  .  _ ^ 

s  -  (0  -  +  ^  ■  f  ^  tan  ) 

-  (« -  iT  f2  -  1  )  -  l«  -  f(M)l 

’  J 


ehsr.  f(M)  1.  given  by  (30.26)  «.d  Table  1.  Consequently  the 
tlsmsnn  invariants  are  v«ry  slsply  expressed  for  the  characteristic 


curves  In  the  physical  plane  as 

l®-f(M)l  -  constant  on  C+  [®  +  fWl  *  constant  or.  C.  (30.36) 

The  function  f(M) ,  ehlch  was  derived  froa  the  trace  of  the 
velocity  veetor  fpllowing  the  characteristics  In  the  hodograph 
plan,  in  (30.26).  Is  isen  to  have  Important  physical  significance, 
md  directly  provides  the  Rlemann  Invariants  for  the  steady  two- 
dlnenilonal  potential  flow.  In  gas  dynamics  f(F)  .  f(M)  Is  refaired 
to  aa  the  Prandtl-Meyer  expansion  function  and,  In  the  form  In  which 
It  la  given  in  Table  I.  It  corresponds  to  the  supersonic  free  expan¬ 
sion  about  a  sharp  comer  as  shown  In  Figure  44  for  the  centered 
siaple  wave  with  a  specific  heat  ratio  of  2.  Since  this  Prandtl 
Meyer  function  la  so  Important,  let  u.  r.-darlve  it  on  another  basis 
chat  will  further  Illustrate  Its  physical  significant-.  From  the 
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fact  that  f  foras  tha  Riaouinn  invariant »  or  from  tha  nature  of  F 
or  F*  alone,  it  ia  evident  that  u  and  w  cannot  be  independent 
of  one  another  on  any  a^ich  ainple  char  ac  ter  i  a  tic*  Conae^uently, 
if  we  write  the  original  potential  equation  (30*13)  in  the 


phyaical  plane  aa 

.2 


and  introduce  w  -  w(u)  ao  that 

Wy  -  UyW*(u),  u^w'(u)  -  -  9jiz 


or  is  . 

dtt  •  « 


yv  ■“ 


7 


-1 


(30.24') 


(  ^  ^  ) 

Thl*  derivation  gives  exactly  the  aamc  result  as  In  (30.24)  and 
verifies  the  fact  that  discontinuities  can  occur  In  the  first 
derivatives  normal  to  a  characteristic  curve.  If  w*  introduce 
(30.32)  Into  (30.24)  v*  obtain  the  equivalent  of  (30.26) 


^(Is )  - 


X 


(30. i7) 


vhlch  again  has  f(P)  as  th*  general  Integral  because  (30.28)  shows 

that  dO  ^  d(q/C*)  _1 _  dF  (30.38) 

F*  ■ 


F^-1 


m 


However,  neither  of  theae  methoda  givea  the  direct  proof  that  f(F) 
providea  the  HLeMrm  invariant*  Thia  fact  «ay  be  proved  directly 
by  the  following  derivation  which  utilizea  the  velocity  component 
0^  along  the  characteriatic,  and  0n  (30*30),  noijul 
to  aa  ahowo  in  Figure  41,  ao  that 
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0  ^  "  q  cos  o(  ,  0n  “  c  ”  q  sin  o< 

d(0^  )  -  (cos  cy  dq  ~  q  sin  (X  d  o<^)  -  c(d9  -  cM 

A  2  -  l/8ln(9+o('),  yU^  -  l/3in(Q-o<) 

since,  from  (30.37),  ^  .  sin  c< 

*  cos 


Then  from  (30.28)  and  (30.32)  we  have 


2 


..2. 


+  1 


-  l/3(q^  ■  <*1  )  “  ^ 


msx 


(30.39) 


(30.40) 


80  that  (30.39)  may  be  written 


(0  -cx)  +  constant  ■ 


-  f3  aln' 


d(gAJ.  .  I  d(gAj 


-■/3 


^(0>  /<lmsx) 


l/l/3(qLx-  ‘•x  >  /W> 


tan 


,r^ 


(30.41) 


Ttils  may  be  finally  written  In  terms  of  (F)  alone  by  noting  from 


(30.39)  that 

li.  -  -  q  -\/f2  , 1,’ 

c  9^  q  sin  cx  tan  oe  i 


consequently  (30.41)  reduces  to 

(«(F)  -  f?  tan'Vfspi  y.  can'^  p  constant]  - 

(«(F)  -  fT  tan-^f  -  I’l  ■  ■  constan'  130.42) 

where  f(F)  Is  the  same  Prandtl-Meyer  function  as  given  In  (30.26) 
and  Table  I.  Therefor,  w.  have  proven  that  the  Rlemann  Invariants 
ere  given  by  (30.36)  and  (30.26).  In  addition  to  the  relation 
betwMn  f  and  P  In  (30.26)  It  Is  aometlmet  convenlant  to  use  on. 


of  the  following 
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f(<v)  -  I  t«n  Of)  +  cx  -  - 

-  f<0  -  ifT '  J  (»•«') 

Th*rii£or«  a  nuMrical  solution  can  now  be  obtalped  for  the 
ganaral  problaa  having  both  families  of  character!  a  tic  a  represent 
curved  non-eieple  waves  by  carrying  on  a  simultaneous  finite 
difference  solution  in  the  physical  (x,ar)  plane  with  (30 .33)1  aid 
in  the  hodograph  (u,w)  plane  by  (30.26),  (30. 3A),  and  (30.36). 

Almost  any  initial  or  boundary  value  data  can  handled  in  this 
manner  as  long  as  the  curve  on  ^ich  the  data  is  given  is  not 
coincident  with  a  characteristic  curve.  The  solution  cannot  be 
obtained  in  the  neighborhood  of  any  portion  of  the  boundary  value 
curve  that  happens  to  be  tangent  to  any  characteristic  curve,  because, 
as  proven  by  (30.20),  the  solution  is  indeterminate  for  boundary 
value  data  given  on  a  characteristic.  It  la  easily  seen  by  this 
finite  difference  method  that  the  dat*  along  a  smooth  non-character¬ 
istic  curve  can  only  determine  the  solution  Inside  the  quadrilateral 
fonned  by  the  characteristic  curves  passing  through  its  end  points 
(Figiire  43)  e.g.,  see  Prelswevk  [1938].  This  well-known  behavior 
of  hyperbolic  type  partial  differential  equations  is  most  directly 
demonstrated  by  writing  them  in  their  normal  or  canncnical  form  by 
transfoming  the  coordinates  to  curvilinear  axes  which  are  the 
characteristic  curves  themselves.  For  example,  Prelswsrk  [1938] 
transforms  tha  aqulvalent  of  (30.23)  onto  the  curvilinear  character¬ 
istic  coordinate  (A,/*-)  systmi  to  obtain 
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£<F^)  »  1/2  (Ay4,),  9-1/2  (Ay*) 

X"  *  ^  ^  (30.43) 

Aic  / 

K<  A  .  A.)  -  <1  -  ^*/2) 

/  . . .  '  •  *7  I  I 

y3’(3-F|)l/2(F2  .^.3/2 


This  normal  or  cannonical  form  Is  so  useful  in  carrying  out  the 
finite  difference  method  of  solution  that  the  values  of  K  have 
tlfo  bees  Incliided  in  Table  !•  It  could  be  used  in  the  following 
type  of  approximation »  as  indicated  in  Figure  43  where  (1,  3)  are 
known  values  and  (2,  4)  are  to  be*  calculated, 


V  -  ~  i/ 

A,  -  A,  '  V 


^  2  ■  ^  1 


{K*  \)  -  X 


(^3 


A 


(30.43*) 

)  -  (^4  ^  2) 

(^4  "  ^ l)  (^2  " 


Consequently  if  the  data  were  given  on  only  one  characteristic 
curve  the  method  would  fail  since  the  values  must  be  known  on  both 
characteristics,  or  on  the  non^ characteristic  curve  s  in  Figure  43 
80  that  one  con  also  wj.lte 

^8  ■  /<8  ■ 


The  ntmarical  method  of  solution  by  finite  differences  following 
(30*43*)  is  known  as  the  ''lattice  point  method"  and  replacas  the 
original  partial  differential  equation  (30.43)  by  a  set  of  linear 
slgtbraic  equations.  The  othei  coosonly  used  semi -graphical  method 
of  solving  hyperbolic  partial  differential  equations  is  callad  the 
network  or  "mesh  method"  and  can  be  Illustrated  by  writing  (30.43) 
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in  thm  form 

A{\.)  -  -  K(  ^;>+  ^X*)  Ak 

^  r  /  (30.43" 

^  ()ju)  -  - 

The  av«rag€  value  at  the  center  of  each  mesh  formed  by  the 
characteristic  network  is  used  for  the  trial  and  error  nuaerlcal 
calculation  of  each  A  increment.  The  increments  are  drawn  tangent 
to  the  characteristic  curves  as  indicated  in  Figure  43.  The 
simultaneous  semi -graph! cal  solution  must  be  carried  out  in  the 
physical  plane  Mb  shown  in  Figure  43  by  using  (30.33)  and  writing 
(30.39)  and  (30*41)  in  finite  difference  form. 

As  a  further  aid  to  nvanerical  and  graphical  solutions  it  is 
useful  to  plot  f(F^)  from  (30.26')  or  Table  1  on  the  hodograph 
(u/c^,  w/c^)  plane  as  shown  in  Figure  43.  The  single  curve  defined 
by  Table  I  may  be  drawn  and  then  rotated  by  equal  increments  of  G » 
or  the  construction  may  be  accomplished  entirely  by  graphical  means 
as  indicated  in  Figure  43  by  rotating  the  small  circle  upon  the 
inner  unit  circle  representing  critical  flow,  while  the  outer  max¬ 
imum  circle  has  a  radius  of  VT  representing  q  (30.29). 

This  geometrical  construction*  yields  f(F^)  since  it  is  an  epicycloid, 
t\B  proven  by  Preiswerk  [1938],  or  Courant  and  Friedrichs  [1948]  page 
262.  All  simple  waves  must  follow  the  characteristic  epicycloid 
in  the  hodograph  plane  because  simple  waves  are  defined  by  (30.37) 
which  has  been  proven  to  have  f(F^)  as  its  integral.  It  ern  be 
shown  that  all  streamlines  corresponding  to  non-simple  waves  must 
lie  within  the  corresponding  characteristic  epicycloids  as  indicated 
in  Figure  43,  since  the  streamline  must  have 
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tan  o< 


(30.44) 


unless  a  finite  disccmtintiity  corresponding  to  a  hydraulic  jump 
(or  shock  wave  in  a  gas)  is  formed. 


Another  useful  aid  in  the  hcdograph  graphical  construction  is 
the  velocity  ellipse  which  is  also  drawn  in  Figure  43-  Wherever  the 
velocity  vector  q  touches  the  curve  of  the  ellipse,  it  will  be  found 
that  the  major  axis  of  the  ellipse  is  in  the  direction  of  the  tan¬ 
gent  to  the  corresponding  characteristic  (either  or  C  )  in  the 
physical  plane  becaus^i,  as  a  consequence  of  (30.29)  and  (30.32),  if 
we  assume  that  9  -  oC  then 

(^f)  -  ) 

fc~)  "  (I  -  8ln^o<)  -  2/2  (F^  -  1) 

\  2  * 
i  §{y  M  -  f2  -  [3  .  2(^)^] 


(u/c.)2 


(w/c*)2j  -  1  j  . 


-  (0»  /c*' 


4  / 


(30.45) 


This  gives  the  velocity  ellipse  shown  in  Figure  43  with  a  major 
axis  of  fY'  and  a  minor  axis  of  unity.  The  major  axis  is  always 
at  the  Mach  angle  o*  with  respect  to  the  velocity  vector  q  because 

we  find  from  (30.28)  and  (30.32)  that  when  © 

^  v2  /  w  V  2  /  v2  A, 


(k)  -ft)  -f^)  -  u-wft)) 
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As  in  th«  previous  csss  of  unsteady  one-di*enslonsl  flow 
over  s  fist  bottoa  we  esn  obtain  very  simple  solutions  for  the 
case  of  sliaple  waves.  In  this  case  thera  is  an  analogy  between 
the  (t.st)  dlagrma  and  the  (x.z)  diagram,  e.g..  see  Courant  and 
Friedrichs  [1948],  As  before  the  simple  wave  corresponds  to  the 

characteristics  of  one  family  becoming  straight  lines,  as 
In  the  exmaples  shown  In  Figure  44,  so  that  (q,  9,  ,  ^  )  *^0  all 

constant  on  the  straight  line  dx/dx  -  constant  In  the  jtoyslcal 
plane.  Therefore  any  given  straight  charactaristlc  line  has 
all  of  Its  properties  determined  by  f(F^)  fto*»  (30.26)  and  each 
of  the  straight  lines  in  the  physical  plane  maps  onto  a  single  point 
of  the  same  single  characteristic  epicycloid  in  the  hodograph  plane. 
The  characteristics  of  the  other  family  remain  curved  in  the  physical 
plana  and  map  In  a  unique  continuous  manner  upon  the  corresponding 
characteristic  epicycloid  arcs  In  the  hodograph  plane.  As  before, 
in  a  simple  wave,  these  curved  characteristics  are  not  required 

for  Si  nuBsrlcal  solution. 

CoMon  exmaples  of  simple  wave  problems  are  shown  In  Figure  44, 
sMd  they  elways  occur  idienever  a  region  of  constant,  uniform  properties 
adjoins  a  region  having  any  variation  in  Its  properties,  the  two 
regions  always  being  Joined  by  a  straight  line  physical  characteristic 
.  constant)  as  long  as  no  finite  discontinuities, corresponding 
to  hydraulic  Jmsps  or  shock  waves,  have  been  formed.  These  finite 
discontinuities  corresp<?ad  to  an  envelope  of  the  straight  character¬ 
istic  lines  that  must  form  whenever  the  boundary  surface  curves  to¬ 
wards  the  oncoming  flow,  resulting  In  a  flow  compression  cr  decrease 
of  velocity  and  increase  In  water  depth  as  Indicated  In  Figure  44. 

The  solution  Is  no  longer  single  valued  at,  or  downstream  of  the 
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MT/elope  so  this  region  must  be  replaced  by  a  hydraulic  Jmp  having 
a  finite  discontinuity  and  non-isentropic  flow. 

If  the  local  flow  velocity  and  water  depth  are  required  only 
on  the  curved  boundary  itself,  then  neither  family  of  characteristics 
have  to  be  determined  (except  as  a  precaution  to  verify  that  no 
finite  discontinuities  have  formed  near  the  boundary  due  to  flow 
compression) *  The  solution  on  the  evurved  boundary  itself  is  given 
directly  from  Table  I  by  simply  measuring  f(F^)  as  the  value 
corresponding  to  (see  Figure  44) 

f[F*(e)l  -  flF*^]  +  0  (30.46) 

If  this  expression  becomes  zero  it  signifies  that  the  supercritical 
flow  has  been  compressed  to  critical  speed  and  a  detached  hydraulic 
Jtaap  can  occur  as  in  Figure  44. 

Vfhenever  disturbance  waves  enter  along  both  families,  either 
due  to  another  boundary  or  by  reflection  fron  a  hydraulic  Jump,  as 
in  Figure  4.2,  then  the  mixed  region  contains  non-simple  waves,  and 
only  a  nxinerical  solution,  similar  to  the  ones  discussed  In  conjunc¬ 
tion  with  (30.43),  can  yield  the  exact  solution*  However,  an 
approxiaiate  solution  for  the  particular  cases  shown  in  Figure  42  can 
be  obtained  by  approximating  Che  curved  cha-scteristlcs  in  the  non- 
simp  le  region  by  means  of  simple  wave,  straight  characteristic  lines. 
The  geoastrical  construction  assumes  that  the  curved  boundary  wall 
of  the  nozzle  can  be  replaced  by  a  series  of  straight  chord  lines 
that  each  have  the  seee  magnitude  of  0  at  every  comer  as 
depicted  in  Figure  45*  At  each  expansion  comer  it  is  assiaed  that 
the  ctntared  simpla  wave  (corresponding  to  a  portion  of  the  complete 
Frandtl*Heyar  expansion,  f)  can  be  epproxioeted  by  a  single  physical 
characteristic  that  is  the  averegc  of  the  actual  expansion  fan  of 
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cb*sr«cterl«tlc«.  Thli  Is  ths(dz/dx)  strslgibit  11ns  thst  is  notMl 
to  tbs  nldpolnt  of  tOis  A  9  splcyclotd  src  roprsssntlng  the  expansion 
angle  ohangs  at  this  corner,  as  shown  In  Figure  45.  Similarly,  the 
cospression  comer  that  turns  Into  the  flow  is  represented  by  the 
single  compression  slstple  wave  that  is  noinal  to  the  midpoint  of  the 
A  9  epicycloid  arc  representing  the  compression  angle  change  at 
this  comer.  It  will  be  shown  that  the  angle  of  this  slngln  av«r*g* 
comprassloo  wave  Is  sctually  the  correct  llsiltlng  value  for  a  weak 
hydraulic  Jiap.  The  geometrical  construction  is  carried  out  in  the 
manner  Indicated  in  Figure  45.  Whenever  a  streamline  crosses  one  of 
these  finite  asiplltude  construction  characteristics  the  flow  is 
Msiaad  to  bend  through  the  A  0  astioclated  with  the  finite  comer 
bend  which  supposedly  produced  this  single  finite  wave.  The 
corresponding  construction  in  the  hodograph  plana  transfars  to  the 
epicycloid  arc  that  is  normal  to  tha  single  finite  wave  in  the 

physical  plana  as  shown  in  Figure  45. 

Also  shown  in  Figure  45  are  the  geometric  construction#  required 

for  the  reflection  of  these  simple  finite  waves  in  the  physical  plane 
from  either  solid  boundaries,  or  constant  water  depth  free  boundaries. 
In  ths  rsflectlon  from  a  solid  boundary  the  original  boundary  slope 
ia  again  attained  by  the  velocity  vector  after  passing  through  ths 
raflactad  wave  which  ben  the  same  strength  for  flow  deflection  as 
tha  original  oncoming  finite  almpla  wave.  In  the  hodograph  plane 
tha  atroMlina  has  gone  from  on#  fassily  of  epicycloids  to  ths  other, 
ending  at  tha  same  value  of  0.  The  completed  aclutlon  for  the  flow 
inside  a  varying  width  channel  having  auparcrltlcal  flow  (F  >1)  la 
prasantad  in  Figure  42.  For  additional  details  and  aids  on  tha 
graphicsl  cooatructlcna  saa  Pralawerk  (193d].  A#  another  axsmple 
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In  Figure  45,  conoiiler  the  reflection  from  e  free  Jet,  hydraulic 
j«np,  or  any  constant  water  depth  free  boundary,  which  must  occur 
In  svich  a  maxmer  that  the  same  water  depth  la  maintained  after 
passing  through  the  reflected  wave  <Ailch  Is  not  only  on  the 
opposite  family  of  epicycloid  arcs,  but  now  must  have  the  negative 
algebraic  strength  of  the  flow  deflection  of  the  original  oncoming 
wave;  consequently  the  value  of  A  9  Is  exactly  doubled  after  passing 
through  the  reflected  wave.  That  Is  to  say,  unlike  the  ordinary 
raflaction  ttao  a  solid  boundary,  the  reflection  from  a  constant  depth 
free  boundary  results  .In  the  opposite  type  of  wave,  an  expansion  wave 

b«coiM8  a  c«i!W«88lon  wave  and  vice  versa. 

In  conclusion  It  must  be  noted  that  this  two  dimensional  steady 

flow  ai^alysls  Is  only  valid  for  a  flat  horizontal  bottom  as  was 
already  shown  by  (29.3)  for  the  linearized  equations.  If  the 
bottom  slope  varies  then  the  Rlcmann  Invariants  do  not  exist,  simple 
waves  do  not  occur,  and  the  numerical  solution  Is  mucl.  more  complicated. 
However,  there  Is  an  even  more  important  criterion  that  must  be 
satisfied  before  ajx  of  the  solutions  given  so  far  can  be  applicable. 
This  is  the  necessary  requirement  that  all  the  perturbation  quantities 
involved  (u  -U,  w,  -h )  must  be  sufficiently  small  so  chat  It  Is  not 
necessary  to  Introduce  the  second  order  terms  from  section  31.  A 
satisfactory  evaluation  of  this  criterion,  at  least  for  F  not  too 
near  unity,  can  be  obtained  by  comparing  the  solutions  of  the  non¬ 
linear  equation  (30.17)  with  the  linearized  equation  (29.23)  or 
(29.32).  As  Is  well-known  In  gas  dynamics,  and  Is  apparent  by 
inspection,  (29.23)  and  (29.32)  are  not  satisfactory  for  F  approaching 
unity  since  additional  terms  must  then  be  retelned.  For  example, 
on  the  boundary  profile  Itself,  (29.3)  for  a  flat  horizontal  bottom 
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iBuat  Incluia  tb«  additional  tam  3F^  (0,^/0)  0^,  which  corraaponda 
to  th«  "transonic  approximation"  of  the  gas  d3rnamic  equation  (with 
a  specific  heat  ratio  of  2),  in  the  limit  as  P  approaches  unity# 

However )  for  the  solution  of  the  steady  flow  everywhere  about  a 

two-dimensional  profile  it  may  be  necessary  to  use 

t(l  -  0^  +  0„1  -  f2[3  ^  0^*  +  2  0^J  (30.47 

alnca  (29.29)  Indtcatas  that  far  from  the  profile  (w/U)  - 

«(0^/U)^l/z  wharaa.  1/^2. 

In  any  case  any  radical  increase  in  the  order  of  magnitude  of  any 
perturbation  term  Isaaediately  indicates  that  the  second  order 
terns  diacuased  in  section  31  must  be  introduced,  since  the  non¬ 
linear  equation  (28# 1)  and  all  the  preceding  results  are  based 
only  on  the  first  order  terms  of  the  shallow-water  theory  a 
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30"^  •  One-iDiMnsioQal,  Steady,  Open  Channel  Hydraulics  and 
the  Hydraulic  Jump. 

The  relations  given  in  eciuations  (30.27),  (30.28)  and  (39.29), 
and  shoim  in  Figure  42,  can  be  used  in  what  is  commonly  knom 
as  the  steady  "one* dimensional"  hydraulics  of  open  channel  flow. 

Here  we  assume  that  even  though  the  channel  width  b(x)  is  varying, 
still  the  valxies  of  q(x)  and  ^  (x)  do  not  depend  upon  z  and  therefore 
do  not  vary  on  any  given  cross-section.  In  conjunction  with  the 
steady  "one-dizaensional"  concept  it  is  necessary  that  w 0  '^O. 
Consequently  the  basic  equations  to  be  used  for  a  flat  horizontal 
bottom  are  given  by  q(x)  -  u(x)  in  (30.27),  (30.28)  and  (30.29), 
and,  in  addition,  by  thii  "one-dimensional"  continuity  equation 
b(x)d(x)i^(x)  -  A(x)u(x)  -  Q  -  constant  (30.48) 

ffhire,  froui  Figure  42,  d!(x)  -  [hQ  (x)]  ■  [A(x)/b(x)]. 

The  validity  of  the  "one-dimensional"  assumption  can  be  consid¬ 
erably  in  error  if  b*(x)  is  large  since  it  is  obvious  that  in  this 
case  w  or  0  cannot  be  small.  However,  the  "one-dijoensional" 
approximation  gives  siurpriaingly  good  nisserical  values  even  in 
supercritical  flow  ^  the  channel  is  well  designed  as  in  Figure  42  so 
as  to  maintain  the  flow  as  uniform  as  possible.  However  in  super¬ 
critical  flow  the  velocity  over  any  cross-section  remains  uniform 
only  near  the  design  Froude  nuad>er  (F).  Freiswerk  [1938]  gives 
the  calculated  and  measured  water  depths  in  a  Laval -type  nozzle 
(the  same  one  as  duplicated  in  Figure  42)  at  varioxis  supercritical 
Froude  niad>era  (F  >  1) .  Hia  results  indicate  that  "one- dimensional" 
hydraulics  giva  a  satisfactory  approximation^  probably  having  an 
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•rror  than  10  par  cant.)  avan  for  critical  or  supercritical  flow* 

This  Mthod  should  ba  aspacially  useful  for  subcritical  flow  sinc<» 
the  isora  exact  nisaarical  solution  is  now  vary  difficult  to  obtain 
bacauaa  the  siapla  aathod  of  characteristics  is  no  longer  applicable. 

The  aost  useful »  and  obviously  the  most  accurate,  application  of 
"ona-dijaansional"  hydraulics  is  to  the  constant  width  rectangular 
cross-section  open  channel  flow.  In  this  application  the  friction 
affect  of  the  vertical  channel  walls  generally  has  a  greater  effect 
on  the  variation  of  q(x,y,z)  than  would  any  of  the  more  exact  terms 
of  the  complete  first  order  shallow  water  equations (28 •  1) ,  which 
have  been  derived  cm  the  assumption  of  negligible  viscosity  effects. 
Consequently  the  "one-dljmensional"  assissption  that  q  -  u(x)  provides 
a  aatisfactozy  approximation  for  the  constant  width  (b),  rectangular 
cross-section,  vertical  wall  channel  having  A(x)  -  bd(x).  Even 
more  important,  thi^.  open  channel  flow  analysis  may  be  fxirther 
generalized,  with  but  little  additional  difficulty,  to  apply  to  a 
bottom  slope  varying  also  with  x.  The  "one '-dimensional"  continuity 
equatiem  (30.48)  then  becomes 

u(x)d(x)  ■  ^  ■  constant  (30.49) 

where  d(x)  is  measured  vertically  from  the  varying  bottom  as  shown 
in  Figure  46.  The  generalization  of  the  Bem^^ulli  eqiiation  (30.27) 
to  include  extraneous  head  losses  (h^)  other  than  those  due  to 
friction,  and  local  variations  in  the  bottom  contour  y(x)  as  shown 
in  Figure  46,  may  be  written  as 

“  (in.t.i:*)  -  [d(x)  +  " 

«  constant  (30.50) 

This  assumes  that  in  steady  flow  the  work  of  gravity,  through  the 
known  average  slope  of  the  flow,  Is  wholly  spent  on  overcoming  the 
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frictional  ravistance. 

Another  relation »  that  is  neces8a’>*y  for  calculating  the  sudden 
additional  hsad  loss  h|^  in  hydraulic  Junps  or  other  discontinuous 
flow  phenomenon,  is  given  by  the  impulse -momentun  relation,  see 
Kelvin  [1886]^  citajiliklgli  41914] ,  oii  .Bakhmeteff  [1932] , 

( SS )  -  !'« O'*"*-* "  - 

2 

•  (1/2  gd^(x)  +  ]  (30, 51) 


which  Is  constant  across  the  hydraulic  Jimp  over  a  flat  horizontal 
bottom  as  shown  in  Figure  48. 

E<]Uation  (30*50)  with  zero  additional  head  loss  (h^^  ■  0), 
gives  the  '*one-dlxnensional*'  solution  for  the  open  channel  flow  that 
has  no  finite  discontinuities  in  the  flow  itself,  and  either  has  tlie 
hydraulic  frictional  resistance  exactly  balanced  by  the  given 
average  slope  for  steady  flow  (so  if  y  -  0  the  surface  slope  is 
parallel  to  the  bottom) ,  or  the  hydraulic  frictional  resistance  can 
be  approximated  by  the  Chezy  fozmula  for  the  case  of  varying  open 
channel  fldw,  see  Bakhmeteff  [1932]  or  Stoker  [1957]*  A.  useful 
concept  for  nearly  all  solutions  is  the  definitici.  of  the  critical 
depth  dyy,  which  corresponds  to  our  previous  definition  of  critical 
flow  velocity  in  (30.20),  that  is,  with  wst  0  9  we  assime  that 


d  -  2/3  h  -  4  .  ,  mol  ,1/3 
★  o  g  I  IJ  J 


(30.52) 


The  last  relation  for  d^  can  be  obtained  either  directly  from 
(30.27),  or  by  substituting  the  expression  -  y(2^/3)gho' 
(30.50)  with  y  and  h^  both  zero.  Also  from  (30.29)  we  have 
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ss 
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2  +  F- 


f2 

»  * 


2/C  ^2  .  2Ei 


2  +  F 


2  » 


2  . 


si  .  (Q/b)^ 
|3 


gd* 


<30.53) 


As  «»  sx^plSg  if  ife  ^pply  SYpiations  (30.50)  and  (30.52)  to 
dstsrmina  ths  flew  ralationa  b«twsen  stations  (1)  and  (29  in 
Flgiira  c46a  w  obtain,  slnca  is  generally  negligible  for  a 
aaooth  variation  in  72  > 

■  [<S)  * 

m  a  satisfactory  approxlaation  from  "one-dixsensional^  hydraulics 
at  least  as  long  as  (72/^^)  s\i£ficlently  sniall.  It  is  interestr 
Ing  to  note  that  here  is  another  resemblance  to  gas  dynacdcs 
behaviour  since 


<(^)>^-  £.rFi<l 


As  another  example,  if  we  consider  the  hydraulic  jimp  shown 
in  Figure.  46b,  now  we  find  that  a  solution  can  only  be  obtained  by 
using  the  impulse-monentm  relation  (30.51),  thereby  proving  that 
the  discontinuous  change  occurring  in  a  hydraulic  junp  must  result 
In  a  head  loss.  If  the  bottom  slope  is  negligible,  as  indicated 
tn  Figure  46b,  then  the  impulse -eioaentiaa  relation  (30.51)  may  be 
written  with  (Q/b)  •  (u^  d^)  -  (u2  d2)  as  first  given  by  Rayleigh [19 14) 

^  7  r,  _2  ^  (uidn)2 


r 


JgdJ  +  diui 

p2_iaz^ 

1  gdl 


.do.  2 


].  * 


'  •  i ^ 

I"  '(ji)  J 


) 


(30.54) 
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or,  if  we  let  the  actual  rise  In  water  level  be  a  -  (.d2~di) 


ax  2,  1/2 


-  tl  +  3/2(|j)  t  1/2 


whsrc 

» + V  -  If  ■  yrr 


(30.55) 


gdj 


-1]  -  ini'll  +  8F^ -  11  (30.56) 


Similarly  (30e54)  can  also  be  solved  for 

^2  ■  ■  'iij'  -  ‘"#<1  *  if '  ■ ) 


(30.57) 


Equations  (30.54)  and  (30.57)  may  be  multiplied  together  to 


yield 


.  sdl+fl.,  .  ^2 


3  ]/i  +i/i  +  srTi  2  2 


2  +  F 


1' 


i/<®,  >®9  (30.58) 

I  1*  Sit 


The  last  inequality  in  (30.58)  is  obtained  from  (30.50)  and  (30.52) 
by  noting  that  in  any  finite  hydraulic  Jixnp  the  head  loss  must  also 
be  finite,,  so  ^  (30.50)  must  be  written  as 


d 


2  X  u? 


h 

L 


1 

A 


/^+  8F^ 


(30.59) 


so  (30.52)  and  (30.59) 
critical  speed  ratio, 


give  the  total  head  ratio,  and  therefore  the 


(30.60) 
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\  vr+  8Fy  - 1)^  +  ^*“1 

/  (/l  +  S^ -  1)^  (2  +  F^) 

Consequently  there  is  no  direct  analogy  between  finite 
hydraulic  lumps  and  gas  dynamic  shock  waves,  as  pointed  otit  by 
Praiswerk  (1938],  since  in  gas  dynamics  the  well-known  Prandtl 

2  4 

relation  for  normal  shock  waves  gives  u^U2  " 

constant  through  the  shock  wave,  e,g,  see  Courant  and  Friedrichs 

[1948]  page  146.  The  equations  are  siz&ilar  only  for  the  limiting 

case  as  the  hydraulic  jump  vanishes  so  ■  F2  ■  Ij^  d2  “ 

and  h^^  -  h^^.  However  this  limiting  process  corresponds  to  the 

isentropic ,  potential  flow  case  %diere  there  is  fui  analogy  for 

small  perturbations  over  a  flat  horizontal  botttmi,  as  previotisly 

discussed.  Also,  as  indicated  by  (30.59)  the  head  loss  and  variation 

in  c  could  be  neglected  until  the  third  order  terms  beccxoe  Impor- 
* 

tant,  so  for  F^  near  unity  the  first  and  second  order  terms  of  the 
hydraulic  Jump  relations  correspond  to  the  gas  dynamic  shock  wave 
relations  having  a  specific  heat  ratio  of  2.  However,  this  is 
identical  to  the  known  fact  that  weak  shock  waves  may  be  considered 
iaentropic  to  the  third  order  of  approximation,  consequently  the 
hydraulic  analogy  to  compressible  gas  dynamics  exists  only  for  small 

perturbations  in  potential  flow. 

There  is  no  direct  analogy  between  the  finite  hydraulic  jump 
and  the  gas  dynamic  shock  wave  because  the  hydraulic  jump  Ixas  a 
head  loss  that  must  be  included  in  the  specific  energy  equation 
(30.50).  This  head  loss  results  in  a  loss  of  Kinetic  energy  that  is 
no  longer  available  as  flow  energy  since  it  is  converted  into  an 


<1 


(30.60) 

continued 
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iRsignificaiit  teiaperatiirc  rise  in  the  water  itself.  In  the  gas 
dynamics  energy  equiation  the  entropy  increase  through  a  shock 
wave  of  course  corresponds  to  a  loss  of  Kinetic  energy  but  this 
is  converted,  through  the  increase  of  the  temperatiare  of  the  gas, 

Into  an  adiabatic  enthalpy  increase  that  maintains  constant  flow 
energy  through  the  shock  wave,  e.g,,  see  Courant  and  Friedrichs 
[1943]  page  125.  The  most  unusual  effect  of  this  loss  in  flow 
energy  (or  hj^)  in  the  hydraulic  jmp  Is  revealed  in  (30,58)  which 
shows  that  the  flow  velocity  downstream  of  a  hydraulic  jump  is 
always  less  than  the  corresponding  gas  dynamics  case,  which  maintains 
constant  so  U]^U2  -  c*  .  For  example,  in  the  gas  dynamic  case 
when  then  u^/c^  (for  ^  -  2),  and  therefore 

U2/c^ 1/Y3\  However  in  a  hydraulic  Jimp  (30,58)  shows  that 
U2/c^  0  when  c>o  (  or 

The  experimental  investigations  by  Bakteoetef f  [ 19321  have  shown 
that  the  hydraulic  jxiaps  in  a  horizontal  rectangular  channel  are  in 
excellent  agreement  with  the  predictions  of  the  "one-dimensional’* 
hydraulic  equations  (30.54)  through  (30.60).  Bakhmcteff  found  that 
depth  increases  as  high  as  IQ(d^)  were  in  excellent  agreement  v/ith 
(30.56).  However,  he  found  that  for  oncoming  Froude  numbers  less 
than  f T,  (i.e.,  Fi<|[3)  the  profile  of  the  normal  hydraulic  Jxznp 
developed  undulations,  and  the  relative  length  of  transition  became 
Indeterminate  because  the  undulating  surface  made  the  region  of 
parallel  flow  increasingly  remote  from  the  start  of  the  wave  front 
as  indicated  in  Figure  46c.  It  is  Interesting  to  note  that  F^  ■ /T' 
con  as  ponds  to  the  maximum  absolute  elevation  that  a  hydraulic  jv*np 
can  reach  with  a  given  (although  there  is  no  limit  to  d2/d^) , 
since  (30.53)  and  (30.56)  may  be  combined  to  give 
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nhich  At:tAin8  ttft  maxiiinii  elevation  of  abovs  the  channel 

bottooi  only  for  «  fT,  at  this  condition  %ie  have 


3/8,  r 

®1 


ii 

20 


2 

'‘l“2  “  10 

This  shows  that  for  all  tha  undulating  hydraulic  j»inps(Fj 


the  change  in  total  head  is  less  than  5  per  cent,  consequently  these 
Jumps  can  be  approximated  by  the  isentropic,  potential  flow  rela¬ 
tions.  This  is  of  great  aid  in  calculating  the  slnnt  or  oblique 
hydraulic  jimps,  as  shown  in  Figure  46d,  since  the  characteristic 
epicycloid  values  of  as  given  in  Table  I,  may  be  used  in 

the  manner  indicated  in  (30.46)  to  approximate  the  change  in 

F,  or  upon  turning  through  sn  angle  0  by  means  of  an  oblique 
hydraulic  jianp  whenever  (d2/d3^) -►  1.  The  comparison  between  the 
val;ie  given  by  fi?i)  in  Table  I  for  compression  to  F2  -  1,  i» 
compared  with  the  exact  values  for  the  corresponding  oblique 
hydraulic  Jixaps  in  Figure  47.  It  is  seen  that  although  the  gas 
dynamic  shock  wave  is  not  a  satisfactory  approximation  for  ?i  fF, 
still  the  isentropic  potential  relation  f(Fj^)  provides  an  excellent 
approxixaation  for  nmch  greater  values  of  F,  since  the  criterion 
for  oblique  hydraulic  Jimps  is  that  the  flow  component  normal 
the  discontinuity  (Fj^  sin  or) 

Th«  exMt  r«l*tion*  for  th«  oblique  hydraulic  Jmp  are  given 
by  Prelawerk  (1938]  and  can  be  obtained  by  almply  adding  the  same 
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v«locity  ccmponent  (w  ■  ■  ^2)  tangent  to  both  faces  of  the 

liy<innillc  jifciup  as  shown  in  Figure  46d*  This  results  In  the 
following  equations  (which  reduce  to  the  preceding  ones  for  a 

normal  hydraulic  Jianp  by  simply  letting  9  “►O  and  ^  Vj2)i 

r-  _  1  t«n  CK  r,  ■  t«n  o<  ]_  /ui/ 

1  2  tan(<x-  6)  [  tan(«  -  S)  J  [  “  sin  oTJ 

fh.]^  fsln  <x  7^  .  ru2/fs<>2\  1  ^ 

2  "  *^1  ( i^/  j  I  sin  («  -  e;  J 


tan  9 


sin  C5f 


sin  (  o(  -  0) 


(30.61) 


The  last  two  equations  in  (30.61)  clearly  show  how  the  oblique 


hydraulic  jump  approaches  the  same  value  as  given  by  the  isentroplc, 
potential  relation  f(F]^)  at  any  value  of  as  long  as  (d2/d3^)-^  1, 


since  they  reduce  to  the  Isentroplc »  potential  characteristic  curve 


given  by  (30.31)  whenever 
IS  previously  mentioned, 


9 


0  and  ( 


(30.61)  shows  that 


1.  As  a  matter  of  fact, 
che  oblique  hydriullc 


jump  angle  (oc)  can  be  approximated  as  in  Flgxire  45  by 


(30.62) 


that  Is,  by  taking  o(  as  defined  by  the  a  rage  line  between  the 
two  characteristics  in  the  physical  plane,  or  the  line  normal  to 
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the  ialdpolnt  of  th€  corr«iponding  ch«r«cterlctlc  epicycloid  ere 
in  the  hodogreph  plane  (see  Figure  45).  The  close  approximation  of  the 
characteristic  aplcycloid  iis  given  by  f  from  (20.26)  and  Table  I,  to 
the  oblique  hydraulic  Jump  relations  (30.61)  is  strikingly 
ili.u5 crated  nhen  they  are  both  plotted  in  the  hodograph  plane  as 
in  Figure  46d  (see  also  Prelswerk  [1938p.  The  trace  of  the 
velocity  vector  for  the  oblique  shock  wave  of  gas  dynamics, 
generally  called  the  “shock  polar**,  is  also  shown  in  Figure  46d 
for  a  specific  heat  ratio  of  2. 

As  long  as  f  from  (30.26)  is  in  close  agreement  \7lth  the  oblique 
hydraulic  Juap  relations  (30.61),  then  the  problems  involving  the 
interaction  and  reflection  of  hydraulic  Jumps  can  be  closely 
approximated  by  the  same  procedure  as  detailed  previously  for  the 
characteristic  epicycloids  involving  compression  waves  in  Figures 
44  and  45.  Whenever  the  required  flow  deflection  (0)  is  greater 
than  that  provided  by  the  epicycloid  passing  through  F^,  as  shown 
in  Figure  46d-  then  subcritical  flow  follows  the  curved  or 
normal  hydraulic  Junp  as  Indicated  by  N  -  1  in  Figure  47.  Similarly 
N  *  2  defines  the  maximum  flew  reflection  angle  (0)  that  can  occur 
vithcrut  ending  in  subcritical  flow  with  a  curved  or  normal  hydraulic 
JuQp.  In  both  cases  two  cur\'es  arc  shown  for  the  oblique  hydraulic 
Jviap,  one  for  the  turning  angle  0  that  will  make  the  flow  critical 
(F2  *  1),  and  the  other  one  Is  the  umxlimia  possible  turning  angle 
max)  oblique  hydraulic  Jump  at  the  given  value  of  F][. 

The  latter  always  produces  subcritical  flow  (F2  1)  as  indicated 

in  Figiire  46d. 

All  of  the  preceding  results  primarily  hold  for  hydraulic 
Jiops  in  rectangular  cross*section  channels  with  a  nearly  horizontal 
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bottom.  BakhiMteff  [1932]  experimentally  shows  the  various  effects 
of  steepening  bottom  slopes.  He  also  generalizes  (30.51)  so  it 
will  apply  to  any  constant  cross-section  shape.  However,  it  must 
be  noted  that  our  equation  (29.3)  shows  conclusively  that  (30.51) 
which  completely  neglects  the  w  velocity  component,  cannot  be 
applicable  to  channel  walls  that  arc  not  nearly  vertical.  Sloping 
sides  on  a  channel  would  increase  the  vertical  velocity  gradients, 
make  a  normal  hydraulic  Jump  impossible,  and  induce  unsteady  vortex 
motions  • 

It  must  also  be  noted  that  all  of  the  preceding  results  are 
valid  only  for  relatively  small  bottom  slopes,  as  Indicated  by  the 
direct  comparison  of  (30.50)  and  (30.51)  with  (28.1)  and  (29.3). 

When  the  flow  Is  rapidly  varying  because  of  large  changes  In  the 
bottom  slope,  then  the  change  in  surface  profile  curvature  Is  so 
pronounced  that  the  pressure  variation  can  no  longer  be  considered 
as  hydrostatic.  For  example,  over  the  spillway  of  a  dam  the  centri¬ 
fugal  force  due  to  the  streamline  curvature  can  actually  exceed  the 
hydrostatic  pressure,  thereby  leading  to  a  pressure  less  than 
atmospheric  with  flow  separation  or  violent  oscillations.  At  present 
spillway  design  is  based  on  semi-emplrlcal  methods  or  model  tests 
since  no  satisfactory  mathematical  analysis  Is  available. 
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31  -  Higher-order  Theories  and  the  Solitary  and  Cnoldal  Waves. 


It  will  now  be  shown  that  many  of  the  preceding  methods  and 
results  based  on  the  shallow-water  approximation  are  valid  only  if 
the  local  variations  in  water  depth  are  not  too  large,  and  the  average 
or  undisturbed  water  depth  is  sufficiently  small.  The  first  require¬ 
ment  Implies  that  the  solutions  of  the  first  order  nonlinear  shallow- 
water  equations  (28.1)  do  not  greatly  differ  (at  least  for  Froude 
nianbers  not  near  unity)  from  the  linearized  solutions  given  by  (29.3) 
or  (29.7).  The  second  requirement  essentially  demands  that  the  depth 
h  be  much  less  than  the  effective  wave  length  A  in  any  application, 
that  is  (h/A)  -  0(1/10),  in  order  to  reduce  the  effects  associated 
with  the  InfinltQsimal-wave  approximation. 

As  already  discussed  the  Infinitesimal -wave  approximation  gives 
a  direct  variation  of  the  fluid  particle  motion  with  the  distance  below 
the  free  surface,  and  also  has  the  propagation  velocity  depend  directly 
upon  the  wave  length,  as  shown  in  section  15.  There  it  was  proved 
that  the  velocity  defined  by 


can  only  be  considered  a  phase  velocity  while  the  actual  rate  of 
propagation  of  energy  is  associated  with  the  group  velocity  defined  by 


UVh/  X 

slnh  ftTTh/A 


(s«.  al.o  L«nb  [1932]  p«ge  381).  Any  such  variation  will  directly 
Interfere  vd.th  the  applicability  of  the  shallow-water  results. 
However,  as  long  as  h/A  <  1/10  it  is  seen  that  the  phase  velocity 
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and  the  group  velocity  are  both  satisfactory  approximations  to  the 
shallow-water  first  order  result  that  c  ■  { ^  and  is  independent  of 
^he  effective  wave  length. 

This  means  that  if  small  scale  model  tests  are  used  to  simulate 
xesults  appropriate  to  the  shallow-water  theory,  then  the  undisturbed 
^aCer  depth  should  be  less  than  1/10  the  principal  model  dimensions. 
Consequently  if  models  less  than  10  cm.  in  effective  dimensions 
are  used,  then  the  test  water  depth  should  be  less  than  1  cm. ,  so  the 
capillary  ripples  produced  by  surface  tension  must  be  consider^id. 

As  shewn  in  sections  15  and  24  the  effect  of  the  surface  tension  (1) 
Is  to  increase  the  phase  velocity  for  the  short  wave  length  capillary 
xlpples  so  that  (31.1)  is  replaced  by 


■1^ 


2  ^  ^)tanh  (31.1*) 

For  ordinary  water  (at  20°C,  T  -  72.8  dynes/cm.,  f  -  0.998  gm/cm^) 


this  gives  the  interesting  result  that  both  the  phase  velocity  and  the 
Sroup  velocity  are  closely  approximated  by  fgh'  for  all  A  >  2  cm. 
i.f  h^l/2  cm.  However  in  any  small  model  tests  the  surface  wave 
patterns  formed  by  the  capillary  ripples  must  be  Ignored  since  they 
3re  short  wave  length  surface  waves  that  are  never  in  accord  with 
the  long  wave  length  shallow  water  theory. 

Except  for  the  section  on  hydraulic  Junps  the  oreceding  shallow- 
water  results  have  all  been  based  entirely  on  (28.1),  the  first 
approximation  to  shallow-water  theory,  and  this  will  now  be  shown 
to  be  limited  to  relatively  small  wave  amplitudes  even  though  the 
complete  nonlinear  equation  (28.1)  be  used,  and  even  though  the 
bottom  surface  be  flat  and  hurlzontal.  The  second  approximation  to 
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8h«l low-water  theory  will  be  shown  to  inmedlately  yield  particular 
solutions  corresponding  to  continuous  permanent  wave  profiles  of 
finite  amplitude  that  can  be  propagated  without  a  change  in  fora  or 
shape  if  viscosity  effects  are  neglected.  These  permanent,  finite 
amplitude  wave  forms  are  the  cnoidal  waves  discovered  by  Korteti^sg 
and  de  Vries  [1895]  which  reduce,  in  the  limiting  case  of  essentially 
infinite  wave  length,  to  the  solitary  wave  of  Russell  [1837,  1844] 
which  was  first  analyzed  theoretically  by  Boussinesq  [1871,  1872] 
and  Rayleigh  [1876], 

The  second  approximation  to  shallow-water  theory  will  show  that 
the  limitation  of  the  nonlinear  first  approximation  to  relatively 
small  amplitudes  is  primarily  due  to  the  fact  that  the  variation  1.: 
the  vertical  velocity  cannot  be  neglected  as  the  wave  amplitude  is 
increased.  This  of  course  invalidates  even  the  rectangular  channel 
hydraulic  analogy  to  compressible  gas  flow,  since,  as  previously 
discussed,  the  principal  assumption  of  the  hydraulic  analogy  is  that 
the  vertical  acceleration  be  negligible. 

The  third  approximation  to  shallow-water  theory  will  then  be 
presented  to  obtain  new  relations  which  will  predict  the  limiting 
heights  of  the  continuous  finite  amplitude  steady  state  wave  forms 
and  give,  for  the  first  time,  the  complete  second  approximation  to 
the  cnoidal  and  solitary  waves.  It  will  be  found  that  the  pressure 
is  no  longer  hydrostatic,  thereby  violating  the  remaining  principal 
assumption  of  the  hydraulic  analogy  and  the  ordinary  cleasical 
shallow-water  theory. 
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31  c/  -  The  First  and  Second  Approximations  to  the  Cnoldal 
and  Solitary  Waves. 


We  will  now  extend  the  perturbation  method  of  Friedjrichs  [1948], 
which  was  used  to  derive  the  nonlinear  first  order  approximation 
(28.1)  to  shallow-water  theory,  to  obtain  the  second  and  higher  orders 
of  approximations  for  the  special  case  of  the  steady  state  propaga¬ 
tion  of  a  wave  independent  of  z  and  t  over  a  flat  horizontal  bottom 


having  y  *  *^oo  **  constant  as  in  Figure  3T* 

First  we  will  show  that  the  only  steady  state  finite  amplitude 

(o) 

solution  of  the  first  order  eqiiation  (28.2)  is  y  “70"  constant 
and  -  u  -  constant.  This  !»  most  easily  proved  by  substituting 

the  solution  of  the  zeroth  order  terms  in  (10.24)  for  steady  water 
flow  over  a  flat  horizontal  bottom,  namely 


^<o)  _  t/o)  -  0,  -  0, 

into  the  first  order  terms  in  (10.27)  to  obtain 

-  Uq  -  constant,  -  0,  -  -  f  g,  -  7^ 


(31.3) 


■  constant 
(31.4) 


since  °  "  Px°^ 

Consequently  the  only  finite  amplitude  first  order  steady  state 

(o) 

golutlon  must  have  x  *  which  would  only  permit  the  hydraulic 
jvnp  as  a  solution  since  •  0  and  u(°>  -  constant  on  each 

side  of  the  finite  discontinuity  defining  the  hydraulic  jump.  This 
is  in  agreement  with  the  well-known  fact  that  the  gas  dynamics  equation 
or  (28.2),  predicts  that  any  finite  amplitude  disturbance  must  form 
a  finite  discontinuity  which  is  a  shock  wave,  or  hydraulic  Jump,  e.g., 
sea  Lamb  [1932]  pages  278,  481.  However,  chi  second  order  approxima¬ 
tion  of  fhallow-water  theory  (10.33)  does  yield  a  permanent  finite 
amplitude,  steady  state  wave  profile  that  does  not  form  a  discontinuity. 
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These  are  called  the  cnoidal  waves,  as  discovered  by  Korteweg  and 
de  Vries  (1895),  and  the  solitary  waves  of  Russell  [1837,  18AA) , 
Bousslnesq  [1871,  1872],  and  Rayleigh  [18761.  In  order  to  obtain  the 
higher  order  approximations  and  limiting  heights  of  these  waves  It  Is 
more  convenient  to  use  exactly  the  same  non-dlmenslonal  variables  as 
introduced  by  Friedrichs  [1948),  and  also  used  by  Keller  [1948], 
iioiTiely , 

2  *  2h2,  o<  -  uj  X,  ^  •  y/h,  H  -  h4o/h 

u(  • '  ,  /t )  -  u(x , y)  /  -1/  (  Of ,  (w  h) 

y(o^)-  /y  (x)/h,  77(o<.^)  -  P(x,y)/fgh  V  (31.5) 

_ 

■Wx)  -  hy<“>  +U.2  h^yd)  +  h5  y(2)  +  .  . 

the  only  difference  In  notation  being  that  x  and  y  are  now  defined 
as  In  Figure  37,  consequently  the  flat  horizontal  bottom  Is  given  by 
y  -  -h^  or^“  -h^^ /h  -  -H,  and  the  expansion  parameter  <5  -  uo^h^ 

Is  used  as  defined  In  (10.23),  with  (31.5)  replacing  (10.21). 

Introducing  the  transfc^rmation  defined  by  (31.5)  Into  (31.4) 
and  the  corresponding  equivalent  of  (10.33)  we  obtain 


-  0  - 

u  (  «,  ^)  -  Uq  - 


constant,  "  V  ^  "  constant  » 


■  0’  4^^  ■  -  0*  4^^  ■  • 

-  0,  V(2)(y^)  .  u^yd)  , 


0  « 

(2) 


-  0.  Yo)  •  -  -  /^^Vof ; 


(31.6) 
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These  expressions  may  be  Integrated  to  obtain 

-H)]~ 


(11.7) 


d«  -  -  (ugf  +  C)  -  -y<l>(or) 


Y 


n“(^)  r 

*  //  w  -  y  f 

o<  O 


V~2(  Xo)  .  />'o  +  »  1 


u. 


u 


o 


/ 


‘cx. 


The  Identities  in  the  last  equation  show  that  the  solution  for  constant 
is  restricted  to  the  unique  value  defined  by 


u(x,y)  -  xij^  -  l(gh(  y,,  +  H)'  -  \f^o  +  )'  - 


constant 


(31.8) 


which  corresponds  to  the  infinitesimal -wave  propagation  velocity 
(28.3)  and  shows  that  the  steady  state  solution  will  be  in  the  neigh¬ 
borhood  of  the  critical  speed  defined  by  a  Froude  number  of  unity. 
However,  u(^)  *  f(^)  now  provides  a  finite  amplitude  steady  state 
solution  that  does  not  form  a  discontinuity , consequently  the 
behavior  of  the  second  order  shallow-water  theory  is  mathematically 
completely  different  from  the  first  order  (28.2)  shallow-water  theory 
or  the  gas  dynaraics  eq*jations.  The  pressure  variation  is  still 
hydrostatic,  since  /7(^)  does  not  depend  upon  and  only  ‘l/’(2) 
ha»  a  direct  dependence  upon  ^  (m  y) » 

Now  in  order  to  continue  the  solution  and  determine  f(^)  we 

o 

must  Introduce  soo^  terms.  By  following  the  same  procedure  as 
used  in  collecting  the  terms  for  (10.33)  we  obtain  for  the 
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particular  case  of  steady  flow  over  a  flat  horizontal  bottom  the 
following  additional  tenriS  that  are  required  for  completing  the 
second  order  solution; 


°  (31.9) 

V<^>(yo)  -  Uo  y^<2)  +  Jl)  y^(l)  .  ^(2)  y(l)^  v'(3)(.H)  .  0 

7r^^^<yo)  - 

These  expressions  were  first  given  by  Keller  [1948]  and  they  may 
be  directly  integrated  to  give  the  following: 


i  (/^+  -[-(^/2)(^2  +  2Hy3,)f 


°  2^  +  H(  Yo  -  U.f 


+  R(o) 

(31.10) 


“off®*  *  («of  +  C)f 


(2) 


cx 


-  1/6  (Y^  +  3H  -  (y„  +  H)  Roc 

Th«  last  equation  for  gives  the  following  expression  for  the  fc 

tern  in  the  surface  profile 


(2) 


1 


“o  ^(.x)  -  I  u^,f2  +  Cf  +(^1/6X  +  3H  -  2H^)  ^ 

-  (  +  H)R  +  Constant 

while  a  sliallar  expression  may  be  obtained  directly  from  77^2' (  Y  ) 


(31.11) 
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6C][iititi.xig  iCs  relfiulon  lii  (j1#9)  and  (31  >10)  so  as  to  obtain 

-  u„7r(2)(/^)  -  -  Ug  1/2  ]  - 

•  Xo 

r_o  (V^  +  2H  y^)f^  -  u^R(4x)  “  ^  4-  constant^ 

Since  (31.11)  and  (31.12)  must  be  identical,  therefore  we 
may  equate  them  and  find  that  f(o^)  must  satisfy  the  ordinary 
differential  equation 


(31.12) 


\zu^. 


^  (£>2^^  +  C, 

/  v^(:/ 


0 


(31.13) 


after  having  introduced  (31.8)  to  eliminate  Equation  (31.13) 

may  be  integrated  to 


(l/3)u®(f„)^  -  -  C(f)^  +(2/3)u®C^(f)  -  constant. 


(31.14) 


Upon 


noting  from  (31.7)  that  f(«)  ■  u^^^(‘^)  and 


^  f 
C  ^ 


1  + 


^  (1) 
/ 


eil 


it  is  evident  that  (31.13)  and  (31.14)  are  the  same  equations  as 
obtained  by  Bousslnesq  [1871,  1872),  Rayleigh  [1876],  Korteweg  and 
de  Vries  [1895],  Lavrent'ev  [1943],  and  Keller  [1948].  The  physical 
significance  of  each  term  in  (31.14)  was  first  pointed  out  by 
Benjamin  and  Llghthill  [1954],  who  derived  (31.14)  in  an  entirely 
different  manner,  starting  with  the  same  series  expansion  of  the 
stream  function  as  Introduced  by  Rayleigh  [1876].  Benjamin  and 
llghthlJ 1  [1954]  use  the  continuity  equation  (30.49),  the  specific 
energy  equation  (30.50),  and  the  specific  momentian  equation  (30.51) 
to  derive  the  equivalent  of  (31.14),  and  then  they  give  a  very  useful 
discuailon  of  the  mathematical  and  physical  behavior  of  its  solutions. 
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The  appropriate  solution  of  (31,13)  for  the  boundary  conditions 
shown  in  Figure  37  is  given  by  the  square  of  the  Jacobian  elliptic 
function  **cn'’  having  the  modulus  0  ^  k  ^1  and  the  real  period 


4K(k) 


/ 

I 


V2 


de 


j  f 


->  2Tr 


1  -  k^stn^e 

■ 

Substituting 

f(or)  -  -B  cn^(Ao^,  k)  (31.15) 

into  (31.13)  we  find  that  (31.15)  is  a  solution  if,  and  only  if, 

0  <  k  ^  1  and 

B  .  A/3  uW  -  - 

Substltuttng  (31.15)  Into  (31.7)  and  (31.5)  w  obtain 


(31.16) 


2k2-l 


'V  M  “  Jc  BUq  (cn^(Ai«x,  k) 

The  boundary  conditions  In  Figure  37c  then  yield 


1  +  0(  i^) 


(o) 


7 


7  (iIj) 


7o  ^ 


(^j 


O  V  ,2 

k 


•  U)  ^h^Bu-  f • 
I  o  O  V 


«  g 
2-1 


T 


) 


a 


-  (wVbu^),  'Vo  ■ 


y  (x)  -  a  cn^  (AUx,  k)  (31.17) 

Then  upon  introducing  (31.8)  and  (31. lo)  into  ^31. 17)  we  obtain 


(31.18) 
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as  the  exact  second  order  shallow-water  theory  solution  for  the 
first  approximation  to  the  cnoldal  waves  of  Korteweg  and  de  Vries 
11895].  The  remaining  terms  In  (31.6)  and  (31.7)  may  similarly 


he  solved  to  give 


where  cn,  sn,  dn  are  the  Jacobian  elliptic  functions  having  the 
argument  (AWx,  k) ,  with  (Au)  x)  defined  by  (31.18).  It  must  be 
noted  that  0  C  k  ^  1  so  k  can  never  become  Identically  zero  for 


two  reasons.  First,  because  for  k  -  0 
cn^(Au)x,  0)  "  cos^  (A^Ox) 

Is  not  a  solution  of  (31.13)  or  (31.14),  and  second,  because  the 
expansions  given  In  (31.18)  and  (31.19)  are  only  valid  for 


a  0  as  k  0. 

The  limiting  case  of  k  -  1  corresponds  to  an  essentially  Infinite 
wave  length  since  K(l)-^^  ,  and  the  cnoldal  wave  solutions  reduce 


P  (xjfJ).  --h(x)  -  y  ^  0  / a 

f8h„  fitr 


a _ 


f  gh 

ui^  -  1  ^  °  ( 

-  (1  +  fc- )  °  fc) 

y  ^ 


) 


5/2 


(31.20) 
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which  provides  the  exact  first  approximation  to  the  solitary  wave. 

All  of  these  solutions,  for  the  cnoldal  wave  and  the  solitary 
wcve,  are  in  axaci:  agreement  with  the  expressions  first  given  by 
Kortcweg  and  de  Vries  [1859],  pages  430-431,  if  one  neglects  the 
terms  of  0(a/h^)  .  It  will  now  be  proved  that  the  terms  of 
0(a/h^  )  must  be  neglected  in  these  first  approximations  becuase 
the  second  approximations  introduce  additional  terms  having  this 
order  of  magnitude. 

We  can  continue  to  the  next  order  of  approximation  by  collecting 
the  remaining  terms  corresponding  to  £  and  adding  some  of  the  £ 


terms  that  are  necessary  in  order  to  complete  the  solution 

r«>  <y„, 


(3)  .  ./(3) 
(3) 


U^  -  'l/ 
P  o< 


U 


(3) 


oc 


u  u 

O  ex' 


(1)  (2) 

u  u 

r)i 


■  - 


(31.21) 


+  1/(2)  (2)  .  0 

r> 


u  v(3)  +  Jl)/(2)  ^  ;7-(3)  +  ^(2)  1/  (2)  .  ^ 

O'  or  ^  /b 

y «> .  u<“  /<”- /P>  >«>  - 


CH 


>( 


-  H)  -  0 

Now  we  can  combine  the  expression  for  in  (31.10)  with 

that  in  (31.21)  to  write 


u 


i  (/b^  +  2H^)R 


CH  Or 


+  S(c<  ) 
(4) 


7 

J 


(31.22) 


Then  the  expreevion  for  1/'  In  (31.21)  yields 
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Ji, 

■  -/  u^^)(c<./3.)dz  -  -  [(/S+  H)  So.  (31.23) 

+  -  20h3/S  2  +  i6„5j  3h^2.2u3)r^^^  , 

t  A') 

The  bouncUiry  condition  defined  by  the  expression  for  in 

(31.21)  thereby  gives  one  relation  for  that  may  be  written  as 

Uo^(3)(^)  -  ^(^>(4)  -  [  «a)^(2)  -T^(2)y(2),.  [u(2)^Uv^3)y(lJ 

-  xr^C  y  )  -  [u(^)/^2)j  .  (u(2)y(l)]  (31. 2A) 


which  may  be  directly  integrated,  upon  substituting  (31.23)  for 
as  ^ 

u^yj3>  ■  [7^''^^Vo)  u<2>(Yo)y^'^J  ■ 

-  constant  (/o  +  5HyJ  -  20H-y2  +  IbH^)!^^  + 

-^(y^3Hy2  -2h3)R^^  +(Vo+H>S<«>-»^^4^vf+2Hy^)f^^  -u2R-iljfy+ 


constant 


[(1/2X  y“  +  2Hy  )f 


o 


R(<>*  )  1  [u  f  +  C] 
o 

.  \/(3)  m  a\f  k.  1  a 


(31.25) 


Another  relation  for  '  may  also  be  obtained  frcxn  the 


other  boundary  condition  defined  by  tie  expression  for  TT 
(31.21),  namely 


(3) 


y„)  +  y^^^  yo) 


(31.26) 


where  JT  itself  may  be  obtained  by  integrating  the  expressions  for 


and  in  (31.21)  to  obtain 

"ex  /i> 


f7^3)(^.A)  -  -  u^u^3)  _  ^(1)  J2)  .  i/2tt/<2)]'  +  constant  (31 . 27) 

Then  substituting  Jf  from  (31.27),  (31.10), 

^^(1)  -  f,  u(2)  from  (31.10),  u(^)  fran  (31.22),  y^^^  -  -(uof  +  C) 

and  from  (31.12)  into  (31.26)  we  obtain  another  relation  for 

y  ,  naoaly, 
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u  yp> 

O  M  ^  ) 


2  u 

+  XI  s  - 

o 


11 

‘i?  (  yo  +  y3  .  8h3  y^)f^ 


2 


u‘ 


12 


■(y2  +  2H/  )R 


O 


72(y2  +  2Hy^)ff^  +u^fR+^(f^) 


+  (u  f  +  C)  u 
o  o 


+  constant 


(31.28) 


These  t»»  expressions  for  (31.25)  and  (31.28),  must  be 

identically  equal,  therefore  since  u  is  defined  by  (31.8), 
we  find  that  the  unknown  function  R  must  satisfy  the  ordinary 
differential  equation 

c 

3“  R  o<o<,  *  (u~  +  3f)R  +  constant 


r  5 


S  a  -  5H2)f^ry 

[  +  (u^  +  H2)f^^ 


-  l/2(  u'^  -  3H^)ff 

o 


CKcX 


+  o 


2  (f«  +  Tvi-  (f)^  , 

*  o  J 


(31.29) 


the  other  unknown  function  S(<x)  having  been  eliminated  since 
u  -  (  y  +  H) . 

O  O 

When  f(cx)  Is  given  by  (31.15),  then  the  solution  of  (31.29)  Is 


2 


n  » 


R(«)  - 


u 


2k2-l 


1  9  I  4, 

1  -  7  T  /cn^(Ao<  ,  k)  + 

o 


(31.30) 


,  1  h2  \  2  3  k^(l  -  k^)  /,  5  h2  \  3  7 

•'>  (2^  -  1)^  ('  ~^j 

and  (31.11)  or  (31.12)  give  the  £.  terin  of  the  wave  profile  as 

.2  ^  /.  9  \  2 


cn  (Ao<,  k) 


••  (uj 


(31.31) 
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Consequently  the  second  approximation  to  the  cnoldal  wave  profile 
is  given  from  the  preceding  and  (31.5)  as 


‘yj  (x)  -  +  co^h^  +  0(£  - 

"  ^  ■  'Vl  I  1  C"^(Adux,  k^  + 


where 

'7l  -  c  -  i^h^  .  (Aio)^  J  (2k^  -  l)(hUg^)^  - 
-  (A«.)^  f  (2k2  -  +  h^)3 


Then  the  boundary  conditions  shown  in  Figure  37c  yield  the 
relations : 


vl:lch  may  be  solved  to  give  the 


second  appraximatlon 
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vherc  now 


3 

The  remaining  t  terms  from  (31.10)  may  then  be  combined 


with  the  £  terms  from  (31.7),  by  mean^  of  (31.5),  to  give 


(31.36) 

For  the  solitary  wave  we  have  k  -  1  and  essentia'*  ^y 


Infinite  wave  length  so  (31.36)  reduces  to 
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’  htr  •(3/4)(h^  )^  sech^(Acox)[l-sech^(AK>x)]  + 


OO 

+  0 


(*J 


(A(a)x)  " 

. 

°fe‘ )  ^ 

u(x«  v)  ^ 


*co 


1/  r~ /r^T  /^r-^  1  ^  0/^^  ) 

(  -  v“oo/  (  o  \^n^/j  [n^  ) 


5/2 


2  2 

■  fe)  5  (2  ^  +  JJL^^8ech^-A«x)-38ech^(A^xW 


/  a 

K- 


V(x.y) 


^  i  fcj-  ^ 

i  fc)  f  ^  J  8ech^(Awx)  + 

■  1  (wt;}  tl  +  +  2  (Hyy««‘=h''*(A<Ox) 

(f^(^)  (^  +  ^  )  8ech2(Au)x)  tanh  (A<ox)  J« 


(31.37) 


tk: 


.2 


7/2 


oo, 


1  fa. 

-7fc 


'1-6  -3  ^*-^8ech^(Aiox^ 


oo  y 


The  celerity  or  propagation  velocity  (c)  of  a  solitary  wave  Is 
defined  by  (31.37)  as  the  constant  uniform  motion  attained  as  x-^00 


u(  og>) 


i 


gh 


■(  gh 


CO 


/l  +a/3(Kf7)  -  ^3/20Xj^)^  +  O(JJ^)  J  (31.38) 


In  Figure  40  equation  (31.38)  Is  shown  to  be  In  better  agreement 
with  recent  experimental  data  than  Is  the  coninonly  used  Bousslnesq 
(1871)  -  Rayleigh  (1876)  propagation  velocity  given  by 
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The  past  success  of  the  Bousstnesq-Riiylelgh  equation,  as 
opposed  to  t)ie  progagatlon  velocities  derived  by  McCowan  [1S91] 
or  Weinstein  [1926],  as  Indicated  in  Figure  48,  is  easily  explained 
when  one  notices  the  close  numerical  agreement  of  the  coefficients 
of  the  Bousslnesq-Rayleigh  equation  with  the  exact  second  approxima¬ 
tion  given  by  (31.38). 

A  comparision  of  the  second  approximations  with  the  first 
approximations  to  the  cnoldal  waves  proves  conclusively  that  only 
t.he  proper  order  of  (a/h^^ )  must  be  retained  for  each  order  of 
approximation.  For  example,  a  comparison  of  (31.18)  with  (31.35) 
shows  that  a  completely  erroneous  second  approximation  would  be 
obtained  by  trying  to  extend  the  first  approximation  to  include  an 
additional  (a/h^)  term.  The  reason  for  this  is  evident  upon  comparing 
the  first  and  second  approximations  for  in  (31.17)  and  (31.34). 

Each  cuccessive  approximation  directly  affects  all  the  coefficients 
of  the  corresponding  (a/h^^)  terms.  Figure  49  shows  the  effect  of 
the  second  approximation  on  a  solitary  wave. 

Of  course  it  must  be  remembered  that  the  expansion  method  of 
Friedrichs  [1948],  that  was  used  to  obtain  all  the  preceding  results, 
is  applicable  only  to  shallow  water,  or  long  wave  length  wave  propaga¬ 
tions.  However,  this  Is  precisely  the  nature  of  the  solitary  wave, 
especially  if  the  amplitude  (a/h^)  is  relatively  small,  since  Its 
effective  wave  length  Is  essentially  infinite.  The  wave  length 
(  X  *“  2K)  is  infinite  by  the  mathematical  definition  given  by  K(l)"^ 
when  k  -  1,  however,  tny  physical  description  would  lead  to  a 
specific  finite  effective  "wave  length"  since  ^ 

a  solitary  wave  of  any  finite  amplitude.  Also  Frledricha  and  Hyers 
[195^]  proved  that  thla  expansion  method  docs  yield  an  existence 
proof  for  the  solitary  wave,  and  thereby  demonstrates  that  it  will 
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at.  least  provide  asymptotic  descriptions  of  the  exact  solution  for 
the  solitary  wave  problem.  The  convergence  of  the  expansion 
Indicated  by  (31.37)  for  the  solitary  wave  seem  reasonable  in  vlev? 
of  their  excellent  agreement  with  experijnental  data.  The  corresponding 
existence  proof  for  cnoldal  waves  (in  the  neighborhood  of  the 
critical  speed  defined  by  a  Froude  number  of  unity)  was  given  by 
Llttman  [1957].  Again  this  justified  the  Friedrichs  expansion 
method,  at  least  as  an  asymptotic  type  of  series  development.  An 
additional  discussion  of  these  existence  proofs  is  given  in  section  33. 


31  ^  -  The  Limiting  Height  and  Velocity  of  Propagation  of  Cnoldal 
and  Solitary  Waves. 

It  is  interesting  to  note  that  the  pressure  is  still  hydro¬ 
static  for  y  ^  0,  but  is  no  longer  hydrostatic  as  the  bottcn 
(y  ■  "b^)  ie  appraoched  with  the  second  approximation.  Similarly 
the  variation  of  the  horizontal  velocity  component  with  depth  below 
the  surface  becomes  important  in  the  second  approximation  only  upon 
approaching  the  flat  horizontal  bottom.  However,  the  vertical 
velocity  component  is  now  seen  to  be  the  principal  variation  that 
violates  the  basic  assianptions  jf  first  order  shallow  water  theory. 
The  first  approximation  given  in  (31.19)  gives  a  monotonic  variation 
in  ^"(y)  that  is  obviously  necessary  from  physical  consideration  in 

order  to  satisfy  the  continuity  equation.  However,  this  monotonic 

3/2 

variation  in  V(y)  is  of  the  higher  order  of  (a/h^)  so  it  can 
be  neglected  in  the  first  order  equations  (26.2)  as  long  as  the 
resulting  local  variations  in  >7  are  sufficiently  small. 
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The  second  approximation  to  the  vortical  v^elocity  component, 
as  given  in  (31.36),  now  shows  that  the  variation  of  'ir(y)  will  no 
longer  be  monotonic  as  (a/h^)  increases.  This  leads  one  to  suspect 
that  thers  is  a  limiting  value  to  (a/ho<,)  for  cnoidal  and  solitary 
waves.  For  example,  (31.36)  shows  that  in  the  neighborhood  of 
the  wave  crest  where  x  0  so  that 

cn^CAwx)  [1  -  (A^x)2]  »  (1  -  J 

then  (y)  actually  has  a  reversal  in  its  direction  if  (a/h^) 
exceeds  the  value  given  by 


81^2 _ 

9k2  H-  2  / 


(31.39) 


for  any  value  of  y 


> 


0. 


This  limiting  value  can  be  substantiated,  at  least  in 
as  k  — >  1,  by  noting  that  (31.33)  has  a  real  solution  for 
only  if 


2 

2k^  -  1 


a 


\ 


o_ 

oo 


<  5  (2k2  -  1)^ 


7  -  37k^(l  -  k2) 


the  limit 
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leading  Co  a  limiting  value  of 


5k^(2k2  -  1) 

7  -  37k2(l-k2) 


(31.40) 


The  most  interesting  application  of  these  results  are  to  che 
solitary  wave  defined  by  k  -  1,  in  which  case  we  find  from  (31.38), 
(31.39)  and  (31. AO)  that  the  limiting  heights  and  the  corresponding 
total  velocity  at  infinity  are  given  by 
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a  \  ^  3 

}  MAX  ^ 


0.7273  >  y 


0.71A3 


i.284  >  1.281 


(31.41) 


Either  of  these  limiting  heights  would  be  satisfactory  for  a  solitary 
wave  since  recent  experimental  investigations  by  Ippen  and  Kulin 
[1955],  Daily  and  Stephan  [1952],  and  Perroud  [1957]  have  shown 
that  under  properly  controlled  conditions  most  solitary  waves  have 
(a/h^)  <  0.7,  the  maximum  recorded  value  being  0.72.  Not  only  are 
the  limiting  values  given  by  (31.39)  or  (31.40)  in  excellent  agree¬ 
ment  with  recent  experimental  data,  but  they  are  consistent  with  the 

order  of  approximation  involved.  The  (8/11)  value  Is  derived  from 

5/2 

the  vertical  velocity  variation  given  to  the  order  of  (a/h^Q) 
by  (31.36),  while  the  (5/7)  value  corresponds  to  the  terms  governed 
by  or  (a/h^)^  in  (31.32). 

Many  attempts  have  been  made  to  determine  the  limiting  height  of 

a  solitary  wave.  However,  nearly  all  of  the  theoretical  calculations 

have  been  based  on  Stokes  [1880^  page  227j  relation  which  assumes 

that  for  the  limiting  heights  of  any  wave  the  wave  crest  must  form 

o 

a  sharp  peak  or  double  point  having  ai  enclosed  angle  of  120  in 
order  to  reduce  the  relative  local  velocity  to  zero  at  the  crest 
Itself,  e.g.^  see  section  33  or  Lamb  [1932],  page  418.  This  120^ 
enclosed  angle  at  the  wave  crest  was  assumed  by  McCowan  [1894], 

Stokes  [1905],  Gwyther  [1900],  Davies  [1952],  Packham  [1952],  Goody 
and  Davies  [1957]  and  Yamada  [1957].  Several  of  these  values  are 
compared  with  experi2z>ental  data,  and  the  theoretical  values  given 
by  (31.36)  and  (31.41),  in  Figure  48.  It  is  seen  that  none  of  these 
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limiting  hnlghtt  for  solltAry  w«ves  are  In  ts  good  an  agreement  with 
the  experimental  data  as  is  (31.41) «  A  reasonable  explanation  of 
the  failviro  of  the  120^  sharp  crest  wave  to  provide  a  satisfactory 
limiting  hel{^t  for  a  solitary  wave  may  be  obtained  by  noting  Chat 


Korteweg  and  m  Vries  [1395]  proved  that  any  flnice  ampllcuck. 
profile  that  did  not  correspond  to  (31.17)  or  (31.20)  would  not 
be  steady  with  respect  to  time.  Conseqiiently,  (31.37)  defines  the 
only  possible  steady  state  solitary  wave,  and  vdien  (a/h^)  >  (8/11) 
the  vertical  velocity  variation  reverses  its  direction  near  the 
crest.  This  probably  leads  to  an  unsteady  wave  erect  that  breaks 
unsynmetrically. 

Equation  (31.38)  shows  that  the  solitary  wave  only  occurs  in 
supexcritical  flow  since  the  Proude  number  corresponding  to  the 
propagation  velocity  is  always  greater  than  unity.  Its  velocity 
of  propagation  is  always  less  than  that  of  the  corresponding  hydraulic 
Jurap  as  cay  be  seen  by  comparing  (30.55)  with  (31.38),  after  expanding 
it  in  powers  of  (a/d^)  -  , 


(30.55*) 


However,  the  cnoidal  wave  c^  occur  in  subcritical  as  well  as  in 
supercritical  flew,  and  as  shown  by  Benjamin  and  LigSithill  [1954], 
the  unciulatlng  flow  in  the  subcritical  region  behind  a  hydraulic 
juap  produced  at  all  Freude  numbers  less  than  VT* may  well  be 
represented  by  tliese  cnoidal  waves.  The  fact  that  cnoidal  waves  can  form 
in  sid>critical  flow  is  easily  shown,  even  in  the  first  approximation, 
by  writing  the  horizontal  velocity  component  from  (31.19)  as 
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»  r  1  .  1  1  for  all  k  ^  1 

^  [^1 

-  M  +  (1  -  “2)  (31. A2) 

Therefore  (31.19)  shows  that  any  definition  or  the  wave  propagation 

velocity  would  be  subcritical  when  1/2.  Stokes  [1847]  hao 

given  two  logical  definitions  of  the  celerity  or  propagation  velocity 

of  pexrmanent  periodic  wave  forms,  and  each  one  would  define  a 

critical  celerity  corresponding  to  a  different  value  of  k,  varying  as 

1/2  ^  <  1,  the  solitary  wave  (k  -  1)  being  always  supercritical 

fox  a  finite  amplitude •  However,  the  existence  proof  for  cnoldal 

waves  by  Uttman  11957]  is  only  valid  for  average  velocities  (defined 

AS  the  velocity  of  the  vertical  plane  that  would  have  zero  average 

flux  across  it)  that  are  near  critical.  An  interesting  physical 

And  mathematical  explanation  of  these  flow  restrictions  is  given  by 

Banjamln  and  Llghthlll  [1954].  The  main  consideration,  as  shown  in 

Figure  55  and  section  35,  is  that  the  finite  ampiiL  oe  periodic  waves 

corresponding  to  k^  ^  0.9  may  be  better  described  by  usirig  ;  iifinl~ 

tesimal  wave  theory.  This  becomes  necessary  because  the  wa’"'-  lengf^h 

2 

of  the  cnoldal  waves  decreases  rapidly  with  k  when  k  is  near  un 
Figure  55  indicates  that  not  only  must  the  wave  length  be  large  cot-  lared 
tc»  the  water  depth.  In  order  to  satisfy  the  shallow  waCfcr  expiu?'  on 
method,  but  also  the  amplitude  of  the  cnoldal  wave  must  becc*ne 
fxtresaely  small  for  values  of  k  ^0.9,  oi  for  subcritical  flow. 

K.ort<nfeg  and  de  Vries  [1095]  have  also  shown  how  negative 
cnoldal  or  solitary  waves  can  be  for*Md  when  the  water  is  very 
ottallow  and  surface  tenalon  (T)  Is  considered.  Their  correct  fir* 
approximation  may  be  written  as 
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“  t  cn^(Aa>x,  k) 


(Awx)  - 


oo 


-2.  Hjtlhml 

^  h  -^t 


(31.43) 


0*9 


nhitr#  the  negative  algsbrelc  vuluc  Is  Assigned  to  the  surfsce 
profile  ^ihenever 


<v^ 

f  y  » 


j  cw.  for  wate: 


(31.44) 


These  negative  waves  have  very  small  amplitude  and  a  very  large 

wave  length,  but  can  create  a  surprising  particle  motion.  It  is 

interesting  to  notice  that  the  deptlt  of  h^  “  1/2  cm,  wtiich,  if  it 

could  be  maintained,  would  eliminate  both  solitary  and  cnoidal  waves, 

is  the  saaM  depth  found  from  (31.1*)  and  (31. 2)  to  give  nearly  the 

sne  value  of  for  both  the  propagation  velocity  and  tho  group 

velocity  of  infinitesimal  WAVds  (also  see  section  15).  Conseciuently 

the  depth  of  1/2  cm  seems  to  be  the  optiman  for  ordinary  water 

(T  -  72.8  dynes/cm.)  whenever  one  uses  small  models  to  stmulatis 

results  appropriate  to  the  first  order  shallow  water  theory  of 

(28.2),  since  this  particular  depth  minimizes  the  effect  of  group 

velocity  and  variation  with  wave  length  for  the  infinitesimal  waves, 

end  minimizes  the  second  order  effect  <h:e  to  the  existence  of  finite 

Y) 

•aplitude  cnoidal  or  solitary  waves.  However,  the  variation  of  ^ 
must  remain  sufficiently  small  since  a  finite  increase  in  ^  above 
h  =  1/2  cm.  could  still  produce  cnoidal  or  solitary  waves.  Also, 
tht  short  wave  length  or  capillary  ripples  that  will  foTyn  must  be 
neglected  in  these  model  tests. 
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2.27 

-  .uo9 

23 

.268 

1.482 

2.34 

-.108 

24 

,2.56 

lc494 

2.41 

-.126 
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38 
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39 
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40 
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41 
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4.44 
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42 
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43 

.078 
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44 
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45 
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46 
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